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EQUIDISTRIBUTION OF THE CONORMAL CYCLE OF RANDOM 

NODAL SETS 

NGUYEN VIET DANG AND GABRIEL RIVIERE 


Abstract. We study the asymptotic properties of the conormal cycle of nodal sets asso¬ 
ciated to a random superposition of eigenfunctions of the Laplacian on a smooth compact 
Riemannian manifold without boundary. In the case where the dimension is odd, we show 
that the expectation of the corresponding current of integration equidistributes on the 
fibers of the cotangent bundle. In the case where the dimension is even, we obtain an 
upper bound of lower order on the expectation. Using recent results of Alesker, we also 
deduce some properties on the asymptotic expectation of any smooth valuation including 
the Euler characteristic of random nodal sets. 


1. Introduction 


We consider {M,g) a smooth (C°°) compact connected oriented Riemannian manifold 
without boundary, and we study the nonzero solutions of 

(1) - AgCy = A^Ca, 

where Ag is the Laplace-Beltrami operator induced by the Riemannian metric g and A > 0. 
Our geometric assumption implies that there exists a non decreasing sequence 


0 — Ar < A2 A A3 A ■ ■ ■ A Aj... —y -i-oo, 


and an orthonormal basis {ej)j>i of L‘^{M) such that 

Vj > 1, —Agtj = AjCj. 

More generally, we are in fact interested in hnite superpositions of such solutions, i.e. 


( 2 ) 


/g7{a;=1[o,a2](-A,)L2(m). 


Recall that, if we denote by iV(A) the dimension of "Ha, one has the following Hormander- 
Weyl asymptotics m 


( 3 ) 


A(A) = 


\0\g{M) 


2-7rtr (f + 1) 


A^ + 0{K 


n— 


where n is the dimension of M and Vol^ is the Riemannian volume induced by g. In the 
present article, we aim at studying elements in "Ha in the high-frequency limit A —)■ -l-cxo. 
There are several natural ways to study these objects: among others, one can for instance 
look at their norms, their nodal sets, their nodal domains, or their quantum limits. In 
the present article, we will address the question of the geometry of the nodal sets. 


V/ G -Ha, Aff:={xeM: f{x) = 0} . 
1 
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Again, there are many natnral questions regarding the geometry of a nodal set A/} or its 
distribution on M - see for instance the book of Han and Lin [SB] or the recent survey of 
Zelditch [6H] on nodal sets of eigenfunctions. For instance, one can try to compute their 
Hausdorff volume, the number of their connected components, or more general geometric 
quantities like their Betti numbers. 

The simplest question seems to be a priori the computation of the volume. In fact, in 
this case, it was conjectured by Yau [66] that there exist two constants 0 < Cg < Cg < +oo 
such that, for every nonzero solution Ca of ([T]), 

CgA < < CgX, 

where is the Riemannian hypersurface volume. In the case of real-analytic metrics, 
this conjecture was proved by Donnelly and Fefferman [27] - see also [H] for the case of 
arbitrary elements in 'H\. In the C°° case, the lower bound was proved earlier in dimension 
2 by Briining [TB] while in dimension n >3, the best known lower bounds are of order 
as recently proved by Golding and Minicozzi [23] and Hezari, Sogge and Zelditch [1T| 162]. 
Concerning the upper bound in the case, Dong, Donnelly and Fefferman obtained an 
upper bound of order A 2 in dimension 2 [2B1I2B] . and in dimension n > 3, Hardt and Simon 
proved an upper bound of order A'^^ for some c > 0 [39]. Thus, even the a priori simpler 
question of estimating the volume of nodal sets remains to our knowledge far from being 
completely understood in a general setting. 

A natural approach is to ask whether the properties expected for any nodal set are valid 
for a “generic” solution of ([T]) or at least for a “generic” element / of "Ha- Ib order to 
formulate a notion of genericity in "Ha, we introduce the following Gaussian probability 
measure on the space "Ha : 


= d^iA{f) := e 


2 



N(A) 

2 


(I/C\... , 


7V(A) 

with 

j=i 


Remark 1.1. We note that, with this convention, the expectation of the norm is 

[ \\f\\‘L^{M)dRA{f) = 1 - 

J-Ha 


In this setting, Berard gave an “average” version [3] of Yau’s conjecture, i.e. one has 

1 ^ + 0(1)). 

More recently, Rudnick and Wigman also estimated the variance in the case where (M, g) is 
a rational torus [56] or the canonical sphere [65] , and Letendre generalized Berard’s result 
to the intersection of several “random” nodal sets [50]. In [67], Zelditch also showed that 
something slightly stronger than Berard’s result holds. More precisely, he proved that, 
given any smooth function u on M, one has 

(4) [ (rWIM/ll.i^>rf^A(/)- \°'|^“„‘* AY / i..(:t)dVol,(M)(l + 0(1)), 

Jha Vo1(§”) A/n + 2 Jm 
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where /*(5o) is the pullback by / of the Dirac distribution, and (.,.) is the duality bracket. 
This result shows that, if we average over the Gaussian measure, the nodal sets become 
equidistributed in M. 

Using this probabilistic approach, we can in fact say much more on the geometry and 
the topology of the nodal sets. In the case of the canonical 2-sphere, Nazarov and Sodin 
gave large deviation estimates for the number of connected components of a random nodal 
set [52] . Lerario and Lundberg also obtained lower and upper bounds for the expectation of 
the number of connected components when one considers more general families of random 
nodal sets [19|. In [5l|, Nicolaescu estimated the number of critical points of a random 
function / in "Ha from which one can also deduce upper bounds on the number of connected 
components. In [Ell E3], Nazarov and Sodin generalized their large deviation result on 
random spherical harmonics and they studied the number of connected components of J\ff 
for general Gaussian random functions. In our context, if we denote by N{f) the number 
of connected components of A//, their result reads 


V(5 > 0, 


lim /xa 

A^ + CxD 


fenA-. 


N{f) 

Jyi 


CLf. 



0 , 


for some non explicit constant % > 0 depending only on (M, g). Note that it does not only 
give the rate of convergence for the expectation but also a large deviation estimate. It is 
also slightly more precise in the sense that the result remains true if we count the number 
of connected components in some rescaled geodesic ball Bg{x, RA~^) for some large R > 0 
- see Theorem 5 in [61] . 

In [35], Gayet and Welschinger proved that, given any closed hypersurface S in R”, the 
probability to find S in the intersection of A// with a geodesic ball Bg{x, RA~^) is uniformly 
bounded from below by an explicit positive constant - see also [51] for earlier results in the 
case of random polynomials. This allows to deduce lower bounds on all the Betti numbers 
of our random nodal sets. In fact, they also proved that, on any smooth compact connected 
Riemannian manifold {M,g) without boundary and for every 0 < i < n — 1, one can find 
explicit constants 0 < Ci{M,g) < Ci{M,g) < -|-oo such that 


c,{M,g)A^< [ k{Aff)dfiA{f)<a{M,g)A\ 

JHa 

where bi is the i-th. Betti number [33l [35]. Recall that ho^Aff) = N{f). Their result is in 
fact more general than what we claim in the sense that it is valid for any elliptic pseudo¬ 
differential operator of positive ordeiQ m. In the case where n is odd, Letendre [50] showed 
the following asymptotic property of the Euler characteristic xi-^f) of ^ random nodal set: 

(5) ) VoWM) + 0(A”-). 

We recall that the Euler characteristic is given by the alternate sum of Betti numbers. 
Finally, in [57], Sarnak and Wigman recently described the universal laws satisfied by the 
topologies of random nodal sets - see also [T9] . 


^In this case, the asymptotics are of order A ™ . 
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These different results show that more is known on the geometry of nodal sets if one 
only aims at probabilistic results. Even if we do not address this kind of questions in the 
present article, we note that probabilistic approaches have been considered for a long time 
in algebraic geometry - see for instance |13] where Kac estimated the number of real zeros 
of a random polynomial as its degree goes to infinity. In “random” algebraic geometry, one 
aims at studying the zero set of random polynomials of several variables instead of random 
superposition of eigenfunctions, and the spectral parameter A is replaced by the degree of 
the polynomials [HI [60l |59l [101 HB EB [50] . 

So far, we have only discussed the mathematical point of view of random Gaussian 
superposition of high-frequency eigenfunctions but we emphasize that they play an impor¬ 
tant role in the physics literature. For instance, in quantum chaos, one is interested in 
understanding the semiclassical properties of a quantum system whose underlying classical 
system enjoys chaotic features such as ergodicity or mixing. In [7], Berry predicted that 
semiclassical eigenmodes of chaotic systems should exhibit the same behaviour as a random 
superposition of waves. This is known as the Berry random wave conjecture. In particu¬ 
lar, any result on random superposition of waves should provide an intuition on what one 
could expect for a chaotic system. In our framework, the main example of chaotic system 
is given by manifolds with negative curvature whose geodesic flows are exponentially mix¬ 
ing. Motivated by this conjecture of Berry, random nodal sets (and nodal domain^ were 
also extensively studied in the physics literature especially in the last fifteen years starting 
from the works of Blum, Gnutzmann and Smilansky [T3] on nodal domains. The case of 
nodal sets was studied by Berry in [8] where he considered (on 2-dimensional domains) 
similar questions as the ones mentioned above on the volume and the curvature of random 
nodal sets. In the 2-dimensional case, Bogomolny, Dubertrand and Schmit conjectured 
that random nodal sets are in fact well described in terms of percolation processes and 
Schramm-Loewner Evolution models [IT] . 

2. Statement of the main results 

In the sequel, we denote by T*M the cotangent bundle of M and an element of T*M is 
denoted by (x;^) for x G M,^ G T*M. From the point of view of “microlocal geometry”, 
it is often more natural to consider the conormal cycle of A/} = {/ = 0} rather than the 
set Aff itself m Ghapter 9]. Recall that the conormal cycle is defined as follows, for every 
/ in ALa, 

(6) N*{{f = 0}) = {(x;0 e T*M s.t. f{x) = 0,^ = tdj for some t ^ 0} C T*M. 

Note that the above defines the conormal only in the set theoretical sense and we will 
later discuss the central issue of orientation. Outside of its singular points, this subset of 
T*M is a Lagrangian conical submanifold and it contains in some sense more informations 
on the geometry of the nodal sets. For that reason, it sounds natural to us to focus 
on the properties of this set. Before giving more explanations on the relevance of the 
conormal cycle in “microlocal geometry”, we mention, as a first motivation, a celebrated 

^This means the subsets of M where / has a fixed sign. 
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index Theorem of Kashiwara m Thm 9.5.3 p. 385]. In the real analytic case, this result 
expresses the Euler characteristic of A/} as a Lagrangian intersection in T*M between the 
conormal cycle of A/} and the graph of dg for a generic function g on M : 

(7) [ Graph dg] n [A^*({/ = 0})] = x (A/}) 

where the intersection is in the oriented sense. 

To our knowledge, the properties of this set have not been studied neither from the 
deterministic point of view, nor from the probabilistic one. In this article, we aim at 
studying the probabilistic properties of the set A^*({/ = 0}) which are somewhat easier to 
consider. For that purpose, we need to introduce the following subset of 

Da := {/ € "Ha ; 3x G A/} such that d^f = 0} . 

By a Sard type argument - see for instance paragraph 2.3 in [50], one can verify that 
= 0 for A large enough. Thus, for /XA-a.e. / in 'Ha, A^*({/ = 0}) is a smooth 
u-dimensional submanifold of 

T*M := {(a:,0 e T*M : ^ ^ 0}. 

In particular, once we choose an orientation, A^*({/ = 0}) can be viewed as a u-current 
[A^*({/ = 0})] in the sense of de Rham [2^158] . It means that, for every smooth compactly 
supported n-form u on T*M, we define 

{[N*{{f = 0})],u:)= [ 

JN*{{f=0}) 

Remark 2.1. We note that A^*({/ = 0}) has two components for / in Da'- 

A^^({/ = 0}) := {(x, G T*M s.t. /(x) = 0 and ^ = td^f for some t G M±\{0}} . 

We will discuss precisely orientability questions in paragraph 13.11 We will in fact use the 
classical conventions of [221 P- 682]. 

Our first main result shows that this defines in fact an function: 

Theorem 2.2. Let {M,g) be a smooth oriented connected compact Riemannian manifold 
without boundary of dimension n. Then, the map 

(8) / G 77a ^ [N*{{f = 0})] G P;(T-M) 

is integrable with respect to the Gaussian measure /xa- Equivalently, for every test form u 
on T*M, the map 

/G77A^([iV*({/ = 0})],a;)GR 

belongs to (Ha, dp- a)- 

In order to state our second result, let us denote by Qg the Riemannian volume form 
on {M,g), by tt : T*M —)■ M the natural projection, and by 7r*Qg the pull-back of the 
Riemannian volume on T*M. Then, our second theorem shows the following universal 
behaviour: 
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Theorem 2.3. Let {M,g) be a smooth oriented connected compact Riemannian manifold 
without boundary of dimension n. Then, one has 

(9) [/V({/ = 0})|df,A(/) = C„ + o(A”-‘), 

with 

Cn = —TT if n is odd, and Cn = 0 otherwise. 

7rVol(§’^-i) •' 

Equivalently, it means that, for every test form uj on T*M, one has 

f m’{{f = 0})]Md^„{f) = Cn(^L=) f X-Ji,Au, + 0(A”-'). 

Jha \Vn + 2J Jt-m 

We note that an important feature of this result is that it shows a very different behaviour 
depending on whether the dimension is odd or not. In the case where n is odd, our theorem 
shows that the submanifold N*{{f = 0}) becomes uniformly equidistributed on the fibers 
TfM while, when n is even, we only obtain an upper bound of lower order. 


2.1. A brief outline of the proof using Berezin integrals. Our asymptotic formula 
can be elegantly derived by using a representation of the conormal cycle as an oscillatory 
integral over odd and even variables in the formalism of Berezin H [371 [251 [63] - see 
section [5] for a brief reminder. Let us now sketch the principle of our derivation. In a 
geodesic coordinate chart, we hrst introduce, for / in "Ha \ L)\, a map 

G{f) : (x, e, t)eT-MxR* ^y = {f{x),tdj - 0 e 

By classical arguments involving wave front sets, we then represent the conormal cycle as 
a current in the sense of De Rham which is obtained by pull-back operation by G{f) : 

[N’if = 0 )] = f G(fr (SS*'(v)dp' Ad-/... A dy"*') . 

JteR* 

Then we use the formalism of Berezin integrals to write the current 6Q~^^{y)dy^ A dy‘^ ... A 
^] 2 der exponential form : 


6i;+\y)dy^ Ady^...A dy'^^^ = ( [ dy^ Ady\..A dy^^^ 

e-'^^^^p-y+^-dy) dpdii. 


(—2i7r) 


(n + ljn + l) 


Thus, the conormal cycle can be expressed as an oscillatory integral over even and odd 
variables : 


(10) \N-(f = Q)] = —f-—l (f 

{-2t7rp+^ JteK* VA("+ii-+i) / 


Then the calculation of the expected conormal cycle reduces to an evaluation of Gaussian 
integrals except that we are in presence of odd and even variables and this exponential 
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representation of the current dramatically simplifies the combinatorics. Integrating over 
the measure and inverting (formally) the integrals yields that we shall compute 


( [ dpdU. 

+1D+1) \JnA J 

Note that this inversion of the order of integration requires rigorous justihcation which are 
provided in subsection 14.31 We also refer to Remark 14.71 for a discussion on this question 
and the so-called Kac-Rice formula. Finally, we compute the Gaussian integral 



( 11 ) 




'Ha 


using Hormander pointwise Weyl’s asymptotics on the spectral projector |12] and super- 
symmetric Gaussian integration [5^. It must be noticed that the Berezin integral was 
also used in the works |3l |12l |29] devoted to the statistics of the critical points of random 
holomorphic sections with applications to counting vacuas in string theory. The general 
philosophy is that Berezin integrals allow to elegantly express determinants as supersym¬ 
metric integrals which hugely simplihes calculations. 


2.2. Applications: smooth valuations on manifolds. We will now explain the inter¬ 
est of the concept of conormal cycle by giving some applications of our equidistribution 
Theorem. First, we recall a “microlocal version” of the generalized Gauss-Bonnet Theorem 
of Ghern [22l equation (20) p. 679-683] later generalized by Fu for subanalytic sets [3T1 
p. 832] (see also Dehnition 4.4.1 and Proposition 4.4.2] for a detailed exposition). 

Theorem 2.4. Let {M,g) be a smooth compact Riemannian manifold of dimension n. 
Then, there exists a smooth, closed, compactly supported n-form 9 on T*M, so that 
It‘m ^ = I 5 Vx G M and for any smooth oriented submanifold of codimension > 1, one has 

(12) x(.’f) = ([lV(A')],») 

where [iV* (X)] is the conormal cycle of X and where x is its Euler characteristic. 

Proof. Our Theorem is a microlocal interpretation of the main result of [22] that we briefly 
recall. Let ttq : UM 1 —)■ M denote the unit tangent bundle over M. Then Ghern constructed 
a (n — l)-form 0 on the unit tangent bundle UM dehned in [221 equation (9) p. 676] 
satisfying the following properties 

• d© = 0 ([221 equation (11) p. 677]), 

• “ ® IS3 P- 51]). 

• for any oriented submanifold X of M, if we denote by VNX C UM the unit normal 
bundle of X in M then we have a generalized Ghern-Gauss-Bonnet identity 

~ 0 = ([221 equation (20) p. 679] and [5^1 Proposition 4. 4.2 p. 52]). 

JvNX 

The form —0 is called geodesic curvature form by Park [55] Dehnition 4.1.1 p. 51]. 
Since M is Riemannian, the metric gives an isomorphism T*M TM which induces an 
isomorphism of cones T*M = T*M \ {0} ~ TM \ {0}. By the above isomorphism and 
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the natural trivial fibration TM\ {0} i—)■ UM whose hber is the group (M>o, x), there is 
an isomorphism of cones p : T*M i—)■ UM x M>o. Choose any smooth function 99 > 0 on 
M>o s.t. dip is compactly supported and dip = 1, such function is easy to construct by 
considering y G C^(M>o),x ^ 0 such that xif)dt = 1 and set ip{x) = ^ xif)dt- Then 
consider the n-form 6 = p* (—0 A dip) and let us check it satishes the claim of the Theorem. 
First the integral of over a hber T*M satishes : Jrp.j^O = ^ ~ 

/u.MxR>o(“®) ^dip = /u,m(-©)/r>o = 1- Then the pairing {[N* {X)],9) satishes 
the identity : ([iV*(X)],0) = {-Q A dip) = /u^xxR>o ^ = 

([ -q)( [ dp) = x{X). □ 

\JVN{ X) /^V^R > 0 _/ 

xiX) 1 

Remark 2.5. In our main theorems, we chose to consider the conormal cycle in T*M instead 
of UM in order to make to make some aspects of the calculation slightly simpler. 

From this result, we deduce an alternative proof of equality (jS]) which gives the mean 
Euler characteristic of random nodal sets : 


Theorem 2.6 (Letendre [SU]). Let {M,g) be a smooth oriented compact Riemannian man¬ 
ifold of odd dimension n. Then, one has 

(13) X m ) V„WM) (^)" + 0(A“-T 

Proof. Let 9 be the n-form of Theorem 12.41 Then, one hnds 

[ xmd^Aif) = f (|iV({/ = 0})l,9)dM/) 

Jn, 


since the map / 1 —)■ ([iV*({/ = 0})],6*) is integrable with respect to the measure djj.\ by 
Theorem 12.21 Hence by Theorem 12.31 we obtain an integral on the cotangent cone : 



X (TV/) dpAif) 


a 


A 


T‘M \\/n-\-2 

(- 1 )^ 9 ') A C, 
A 


xeM \JT‘M 


MHg A9 + 
A 


(- 1 ) 


a 


x&M V x/n -\- 2 


n + 1 


2(-lf(-l)-i 

7rVol(S’^ 


yjn -\- 2 

Qg + 0{A^-^ 

A 


+ 0{A 


n—l\ 




where we used a Fubini Theorem on the cotangent cone T*M and we integrated on the 

7 ^1 I 1 TTi 1 

hbers of T*M hrst. Note that since n is odd, (—!)"'(—1)^“ = (—1)^“ = (—1)^“ which 
explains the constant found in the statement of our Theorem. □ 
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In [50], Letendre computed the Euler characteristic by expressing it as the integral over 
{/ = 0} of a certain curvature form depending on f as for instance in [2T1 Eq. (9)]. Here, 
instead of this approach, we use Theorem 12.41 which gives the Euler characteristic as the 
integral of a fixed curvature form over N*{{f = 0}). In some sense, this point of view 
allows us to extend Letendre’s Theorem into an equidistribution result in the same way as 
Zelditch’s equidistribution result (|3|) generalized Berard’s result on the volume. 

Chern’s formula for the Euler characteritic of a smooth submanifold can in fact be 
understood in the more general framework of the theory of smooth valuations recently 
developped by Alesker [2] and that we will now briefly review. We consider X a smooth 
oriented manifold (not necessarily endowed with a Riemannian structure) and we denote 
by V{X) the set of all compact submanifolds of M with corners. In the terminology of [2], 
we say that a map 

f : V{X) ^ C 

is a smooth valuation if it is a finitely additive functional, and if it satisfies certain continuity 
properties. For simplicity of exposition, we remain vague on these two notions which need 
to be defined carefully on V{X) - see Part II of [T] for details. As was already mentionned, 
smooth valuations on manifolds generalize classical concepts from integral geometry such 
as volumes, Euler characteristic and mixed volumes. One of the remarkable property of 
these valuations is that they can be represented as follows : for every smooth valuation 0 , 
there exists a smooth differential n-form on T*X such that, for any P in V{X), one has 

<P{P) = [ 

Jn*{p) 

where N*{P) is the conormal cycle of P [2]. We emphasize that the converse statement 
is also true, and that the n-form is a priori non unique. In the example of the Euler 
characteristic y. Theorem 12.41 gives such a n-form as soon as we have hxed some arbitrary 
Riemannian metric on X. 

Thanks to this interpretation of smooth valuations as integrals over conormal cycles. 
Theorems 12.21 and 12.31 imply the following corollary on smooth valuations: 

Corollary 2.7. Let {M,g) be a smooth connected compact Riemannian manifold without 
boundary of dimension n. Let f be a smooth valuation in the sense of |1] such that up can 
be chosen compactly supported in T*M. Then, the map 

f ^ fiAff) 

belongs to L^{'H\,dfi/fj. Moreover, one has 

(14) [ 4,(Arf)dl‘A(f) = cJ^L=) f 7 r-n,Au,, + 0(A"-‘), 

J-Ha \Vn + 2J Jt*m 

where 


7rVol(§"-i) 


if n is odd, and = 0 otherwise. 
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2.3. Further questions. Let us now mention natural questions that can be asked about 
the conormal cycles attached to nodal sets. 

• In the spirit of [50], it would be natural to consider intersections of nodal sets and 
compute the expectation of their conormal cycle. This could probably be directly 
obtained with the methods of the article but would require some slightly more 
complicated combinatorial arguments. For simplicity of exposition, we decided to 
consider only one nodal set. 

• Regarding Theorem 12.31 it would be natural to understand the asymptotics in the 
case where n is even. At least in the case where the set of closed geodesics is of 
zero Liouville measure, we believe that one should be able to compute the term 
of order in the asymptotics using the fact that in this geometric framework, 
the remainders in Weyl’s asymptotics are of order o(A”“^). Understanding the 
geometric meaning of the lower order term (and also checking that it is a nonzero 
term) is probably more subtle than computing the leading term we obtain here. 

• In this article, we only computed the expectation of the conormal cycle and it would 
be of course natural to look for variance or large deviations estimates in the spirit 
of [811561I5211651153]. 

• Finally, it would be natural to understand what happens in the deterministic case. 
Can one obtain at least some upper bounds on the conormal cycle of a deterministic 
nodal set? If yes, what are the rates, and do they also depend on the dimension? 
This kind of questions would require completely different techniques like the ones 
in [35] . 

2.4. Organization of the article. Section [5] gathers preliminary results that will be used 
in our proof. More precisely, our first task in the present article is to express the conormal 
cycle [A^*(/ = 0)] explicitely in terms of Dirac distributions and differential forms. This 
kind of representation appears in the book of Schwartz [58] and it is in some sense slightly 
more adapted to our problem than the geometrically appealing definitions appearing in 
the works of Kashiwara-Schapira, Fu, Alesker, Bernig ([G] HTl HG] [31]). Thus, we give in 
section [5] an explicit integral formula for the conormal cycle. Moreover, we also introduce 
in this section other relevant tools needed for the proof of our main Theorem. 

In section jT] which forms the core of our paper, we give the complete proofs of Theo¬ 
rems 12.21 and 12.31 using microlocal analysis and combinatorics. 

Then, in section (5] building on tools from quantum held theory, we give a new integral 
formula which expresses the conormal cycle of a nodal set as an oscillatory integral in 
bosonic (even) and fermionic (odd) variables. Our formula makes use of the so called 
Berezin integral [UlS^ESllGS] which was already mentioned above. This formula is inspired 
from a very general formula expressing integration currents as oscillatory integrals on odd 
and even variables used by A. Losev et ah in their works on instantonic quantum held 
theories [301 El]. We give a leisurely introduction to the necessary tools to grasp the 
meaning of our formula. Using this very simple integral formula representing [A^*(/ = 
0 )], we give a fast derivation of the leading term of the asymptotics of = 

0})]d/iA(/) which becomes a simple exercise of Gaussian integration with respect to odd 
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and even variables and completely avoids the heavy combinatorics of the second part. The 
introdnction of even (fermionic) variables allows us to avoid complicated combinatorics 
involving sums over partitions at the expense of a little bit of abstraction. The calculations 
are similar to those encountered in quantum field theory on curved Riemannian manifolds. 
This somewhat alternative derivation of the main asymptotic formula is presented in this 
last part in a purely formal way and it relies on several inversions of integrals and a priori 
estimates on remainder terms which could in fact be justified by using the arguments of 
the second part of the article. 

The article also contains 4 appendices. Appendix recalls some classical results on 
the derivatives of the spectral projector of the Laplacian. Appendix |B] justifies the inver¬ 
sion formulas used in the article via tools from microlocal analysis. In appendix [21 we 
briefly discuss a question of independent interest. More precisely, building on the proof of 
lemma IB. 31 we give a functional analytic proof of Whitney embedding theorem. Finally, 
appendix [D] gives a brief account on the properties of currents needed for our purpose and 
carefully justify the meaning of our formula using tools from microlocal analysis. 
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3. Preliminary results 


In this section, we give the three main ingredients for the proof of our main Theorem. 
The first ingredient is Proposition 13.71 which gives a representation formula for the conormal 
cycle in terms of pull-back of Dirac distributions. As we shall see in this proposition, the 
conormal cycle can be thought as a local functional on 2-jets of functions. Thus, our second 
ingredient (Prop. 13.9p describes the pushforward of the Gaussian measure /ia on "Ha by 
the map 


Jliy) : / e 77 a 



G 


where we identify with the hbers of the 2-jet bundle. The pushforward measure 

is in fact defined by its characteristic function which is expressed in terms of the 
kernel of the spectral projector defining our Gaussian measure Then the third ingre¬ 
dient (Prop. I3.10p is a kind of “Fubini statement” which allows to intertwine the order of 
integration of the variables (?/, rj, t, f) involved in our problem. 
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Remark 3.1. For any smooth oriented manifold X of dimension m, we will denote by 
{Vd{X))Q<d<m the smooth compactly supported d-forms on X and by {V'^{X))Q<d<m its 
dual, i.e. the set of currents of dimension d on X. We refer to appendix [D] for a brief 
account on the theory of currents which is necessary in our proof. 

3.1. Representation of the conormal cycle. Our hrst goal is to represent the conormal 
cycle iV*({/ = 0}) in terms of Dirac distributions. For that purpose, we hrst need to recall 
the dehnition of push-forward in the context of currents. Let A : X i—>■ D be a proper 
smooth map between two smooth oriented manifolds X and D, and let T be a current in 
D^(X), where d < dimD. Then the push-forward A*(T) is dehned by duality from the 
pull-back : 

Vcn e V,{Y), (A,T, u:)y = (T, A*n;)x. 

An important fact concerning the push-forward operation is that an orientation of T 
induces canonically an orientation of A*T. 

3.1.1. Abstract definition of the current of integration. Let us start by recalling the dehni¬ 
tion of the conormal cycle of {/ = 0} : 

Definition 3.2. Let M be a smooth connected compact oriented manifold of dimension 
n without boundary. Let / G if d^f never vanishes on {/ = 0} then we dehne 

the conormal cycle [X*({/ = 0})] as the integration current on the conical Lagrangian 
submanifold {(x;^) s.t. /(x) = 0 and ^ = td^f for some f G M \ {0}} in V'^{T*M). 

A subtle fact we must add to our dehnition concerns the orientation of the conormal 
cycle. First, note that the conormal cycle contains two components 

N^{{f = 0}) = {(x;.^) s.t. /(x) = 0 and ^ = tdxf for some t G M± \ {0}} 

since it does not meet the zero section of the cotangent bundle. Hence an orientation of 
[iV*({/ = 0})] consists of a choice of an orientation for each component separately. In the 
sequel, we denote by M* the set M \ {0}. Let us give another description of [X*({/ = 0})] 
in terms of Lagrange immersion which will be useful in the sequel since we will interpret 
the conormal cycle [iV*({/ = 0})] as the push-forward of a cone S' x M* and the orientation 
of [iV*({/ = 0})] will be induced by an orientation of S x M* : 

Lemma 3.3. Let f G C°°{M) where df never vanishes on {/ = 0}. Set S to be the 
abstract manifold {/ = 0} and i : S ^ M the canonical immersion of S in M. Consider 
the cone S x M*. Then, the conormal cycle N*{{f = 0}) is the image of the immersion : 

(15) A : (x;f) G S X M* i —> {i{x);tdi(x)f) ^ T*M, 
and, at the level of currents, we find that 

(16) [N*{{f = 0})] = [A(S X M*)] = AJS x M*], 

where we orient the current [S x M*] by any differential form a G f2”(S x R*) such that 
df A a|5xR>o = klg f\ dt|5xR>o df A Q;|5xR<o = —^g A dt|5xR<o- 
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Remark 3.4. The cone S x 
should orient separately. 


contains two components S x M>o and S x R<o that we 


Proof. By dehnition, [At*({/ = 0})] is just the integration current on N*{{f = 0}) which 
is equal to X{S x R*). Therefore, for any test form u in 'Dn{T*M) : 




bJ = 


bJ 


'N*({/=0}) 


'A(S’xK*) 


= / X*u= {[S xR*],X*u) = {X,[S xR*],u) 
J{SxR*) 

by dehnition of the push-forward which yields the claim. 


□ 


3.1.2. Expression in local coordinates. Our goal in this paragraph is to show that the 
conormal cycle [N*{{f = 0})] can be written (in local coordinates) as the current 


(17) 




(/. /\ d (is,,/ - , 


where A is any positive number - see Prop. ID. II in appendix [D] for the meaning of the 
convenient notation /A — ^). We need to do this progressively and the hrst 

step is to dehne some integration current [2 ’j a] on R* x T*M whose push-forward along 
the hrst factor R* yields the conormal cycle [Ai*({/ = 0})]. For that purpose, we set 


t 


Zf,A ■= (t; X, 0 e IR* X T’M : f{x) = 0, and ^ = -4/ 


(18) 

The following lemma holds: 

Lemma 3.5. Let f G C°°{M) where df never vanishes on {/ = 0}. Let [Zf^\\ be the 
integration current on the set Zf^A defined by equation l[IE) and n the projection {t;x,f) G 
R* X T’M I —(x,0 e T*M. Then 

(19) lAf({/= 0})| = 

Proof. Let bj{x,^]dx,d^) be an element in V^{T*M). One has : 


(4[4,a]i^) = ([%A],T*a;) = 




u{x,f-,dx,df) 


by variable change. 

Hence {TTfiZf^A\,uj) = a;(x, t4/; dx, d(f4/)) = J^^^fiX*u){t,x;dt,dx) since 

{X*u){t, x; dt, dx) = a;(x, t4/; dx, (i(t4/)), hnally = A*[S'xR*](a;) = [N*{{f = 

0})](a;) by lemma 13731 □ 

The next lemma aims to represent the current [2 ^/,a] by Dirac distributions in a local 
coordinate chart. Let {U, fi) be a sufficiently small coordinate chart centered at a point xq 
in M. It induces a coordinate chart on T*M as follows: 


$ := r> : T*U ^ fi{U) x R", (x, 0 ^ {y, v) = (0(x), 
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We also set 

f*U := {{y,r]) G 0(f/) x M” : ^ 0}. 

In the following of the article, we will use these conventions. In particular, (x, will always 
denote a point in T*M and (y, rj) its image in local coordinates. Using these conventions, 
one has 


Lemma 3.6. Let f G C°°{M) where df never vanishes on {/ = 0}. Then, for any system 
of local coordinates {y; rj) on T*U where U d M is some bounded open subset, one finds 
that, in these local coordinates, 

(20) (/(»). df A ^ d(jd,J - r,.) 


Proof. The submanifold Zf \ in M* x T*U is dehned in the local coordinates {y,rj) by the 
system of (n + 1) equations : 


f{y) = 0 

t df 

Xay-oi = 

The collection of one forms [dyf^dtji — d j is linearly independent along 

as 


dyf P /\{dr]i- d 


2=1 


A dy^ 


= dyf A drji A ... A drjn + lower order terms in drj, 


where dyf A dpi A ■■■ A drjn does not vanish on since df does not vanish on {/ = 0}. 
Then, by Corollary ID. 41 the current [Zj a] is represented in local coordinates {y,r]) by the 
formula : 

(21) [%a]|r*xT*u = jdyf - dyf ^/\d (J^dlyj - rji 

□ 



Then combining the above Lemma to Lemma [3.51 which claimed the identity [A^*({/ = 
0 })] = 7r^[Zf^A], we conclude this hrst part by the first key ingredient of our proof: 

Proposition 3.7. Let f G C°°{M) where df never vanishes on {/ = 0}. Then, for any 
system of local coordinates {y,ri) on T*U where U <Z M is some bounded open subset, 
the restriction of the conormal cycle [A^*({/ = 0})] on T*U can be described in the local 
coordinates [y, rj) by the integral formula : 

(22) |iV({/ = 0})1 |t.i, = (/(a), jdj - dj A Ad[tf,.f - >?.) 

where we use the convention + Jq 
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Our convention for the integral = /o^°° + /o °° reflects our choice of orientation for 
the conormal cycle. We note that we use in fact the same convention as in [2^ p. 682], 


3.2. Pushforward of the Gaussian measure on 2— jets over U. Thanks to fl22|) . the 

conormal cycle of {/ = 0} at a point x eU can be thought of as a functional of the 2-jet 
of the function f at x E U. As we aim at computing its expectation with respect to the 
variable /, the pushforward of the measure pA on the 2-jet bundle will naturally appear in 
our problem. The pushforward of this measure will then “live” on the space M” +*^+1 since 
we are working in local charts. From this point on, we will in fact make use of geodesic 
normal coordinate charts for which one has proper asymptotics for the spectral projector 
associated to the Gaussian measure /ia - see appendix 


Remark 3.8. By geodesic normal coordinate chart, we mean the following. Consider the ex¬ 
ponential map Exp induced by the Riemannian metric on M. There exists a neighborhood 
U oi M X {0} in TM such that Exp induces a smooth diffeomorphism {x, v) i— )■ {x, exp,,,(u)) 
from U onto a small neighborhood of the diagonal in M x M. Then, we £x a point Xq in 
M, and we consider the preimage of under this diffeomorphism. This defines a local 

chart around xq into which can be identihed with R” via a measure preserving linear 

isomorphism. 

We £x [y, rj, t) in T*U x R* and m > 1. We dehne the following map from "Ha fo R"'^+"'+b 



For the coordinates in we set 

F = (Fq; Fi, . . . , Fn] Fi^i, . . . Fri,l, Fi^2, • • • Fn,2i • • • , F 2 ,„, . . . , Fn,n)- 

The next Proposition describes the Fourier transform of the pushforward of the Gaussian 
measure /xa under the map J\{y)- It follows directly from the results on the asymptotics 
of the derivatives of the spectral projector recalled in appendix (more precisely Corol¬ 
lary |A2]). 


Proposition 3.9. The pushforward of the Gaussian measure under the map J‘f{y) 
induces a Gaussian measure h\{F) on R”' +"'+^; 

(23) (.F), 


whose Fourier transform has the following expression: 

(24) f ^ ^-v^^{TAA(y)T) ^ 

where the matrix Aa(|/) is defined as follows : 


A^{y) := Ao + C?(A"') 


Af{y) Af{y)) 


1 
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with the constant in the remainder uniformly hounded in terms of y and A, and 


Aq : — 


^ 1,1 ^ 1,2 \ 
^ 2,1 ^ 2,2 / ’ 


where the matrices Aij do not depend on y and are defined as follows. 

Ai^i is an (n + 1 ) x (n + 1 ) matrix, Ai ^2 is an (n+ 1 ) x matrix, A 2 ^i is an x (n + 1 ) 
matrix, and ^ 2,2 is an x matrix. Their expressions are given by 


Ai,i := 


^ 1 4 _ 4 ^ _ 

Onxl J’ ^2,1 -A,2- ^^2 


5 ^n'^xn j t 


and 


^2,2 ) 


(n + 2)(n + 4) 
with Bg^s' an n X n matrix which is equal to 

Bg^g . Id 2 d%ag{^{^6j^fj\.cj<n) if s s , 


and 

Bs,s' Bgi g . (di^gdj^g/ -|- if s ^ s . 

3.3. Intertwining the orders of integration. In our proof, we will need to integrate 
functions involving the variables {y,p,t) in T*U x M* but also / in "Ha. The fact that we 
can “intertwine the order of integration” plays a central role in the following - see e.g. the 
proof of Proposition 14.41 Recall from paragraph 2.3 in [50] that 

D\ := {/ G TLk ; 3i/ G /~^(0) such that dyf = O} . 

is of /lA-measure 0 for A large enough. One can also verify that Oa = 'H\\D\ is an open 
subset of "Ha- For every / in Ba, we dehne the following map on TLa x 4>{U) x R” x R*: 

Ga : {f,y,r],t) G ?^aX 0(17) xR”xR* 1 —^ (^f{y)-,jdyj - r/i,..., jdy^f - G RxR”, 


Then, one has 

Proposition 3.10. Let 4/(/, y, 77 , t) be a test function in 'P(Oa x T*U x R*), then 


(25) (l,(GA(/)-(C‘).«>T.„«r)„_^ = 

Proof. See section [B] of the appendix. □ 
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4. Proof of Theorems 12.21 and 12.31 


Recall that we aim at computing the expectation of the current of integration on the 
submanifold N*{{f = 0}). According to Proposition 13.71 which gives an integral repre¬ 
sentation of the CO normal cycle, one knows that, for every / in and for every u in 
Vn{T*U), one has 


777 / (^{y,dy,r],dr]) = (GA{f)*{6Q^^),df A /\di-dyj -r]j] Aui{y,ri,dy,dr]) 


T*UxR* 


where (•, is the duality bracket in T>'{T*U x M*) x T>{T*U x M*). Recall that we 

chose a particular orientation for M* in paragraph 13.11 

The main result of this section is the following Theorem from which Theorems I2.2l and l2.3l 
follow by partition of unity: 

Theorem 4.1. ITe use the above notations. For any element 


in 'Dn{T*U), let 


UJ ■.= ^ Ua,p{y:V)dy'^/\dr]p 

\a\ + \ 0 \=n 


|a;|| = sup ^ \u:a^ 0 {y,vi)\. 

o‘,l3,iy,v)e:T*U 


For every compact set /C C T*U, there exists a constant Cjc such that for all test form u 
supported in /C 


'Ha 


'y*({/=o}) 


u{y,dy,r],dr]) 


dfiA ^ Cic\\uj\ 


and 


= a 


Ha \aN*{{f= 0 }) 

A ^ ” 


\/n 2 J Jf*u 

with Cn = 0 if n = 0 {mod 2 ) and 


if n = 1 {mod 2 ). 


uj{y,dy,r],dr])j d^A{f) 

dy^ A dy"^ ... A dy'^ A u){y, dy, y, dy) + (P(A”“^), 


a := 


TT Vol(§"'“^) 


Remark 4.2. It follows from the hrst statement of Theorem 14. II that for every test form cn, 
the map 




/7V*({/=0}) 


uj{y,dy,y,dy) 
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belongs to For any integral current T, let M(T) denote the mass of the 

current T [3S]- Our result is in fact slightly stronger than what we claimed in Theorem 12.21 
as it also means that, for any smooth compactly supported cut-off function ip G V^{T*M), 



M{p[N*{{f = 0})])<+oo. 


In other words, the mass of the cut-off conormal cycle (p[N* ({/ = 0})] is a function of 
feUA. 


The purpose of this section is to prove this Theorem. We emphasize that there are two 
parts in this statement. On the one hand, we have to show the integrability property and 
on the other hand, we have to compute the precise value of the expectation. The hrst part 
of the statement requires a delicate analysis which is carried out in the hrst part of this 
section. After that step, we can combine this hrst part to some combinatorial arguments 
in order to obtain the second part of the Theorem. 

More precisely, the proof is organized as follows. First, in paragraph 14.11 we write a 
formal asymptotic expansion in powers of of 




1 


W*({/=o}) 


u{y,dy,r],dr]). 


Paragraph 14.31 is the more delicate part where we prove the hrst part of Theorem 14.11 
concerning the integrability of the conormal cycle as a function of / G "Ha with respect to 
the measure dpA- After that, in paragraph 14.41 we compute explicitely the expression of 
the leading term and in the remaining paragraphs, we compute the combinatorial constant 
appearing in the leading term. 


4.1. Preliminary simplification. Before getting into the details of the proof of Theo¬ 
rem STl we start by making a few reduction that will make the calculation slightly simpler: 

Lemma 4.3. Let 

ui := ^ uja,yiy, V)dy°‘ A dpp 

\a\ + \y\=n 

he an element in Vn{T*U). With the above conventions, one has, for every A and for every 

f in 

^ [ ^{y,dy,V,dr]) 

A" ijv*({/=o}) 

n-l 

k=0 


for every 0 < k < n — 1, Pk{R, S,T) is a polynomial which does not depend on A, 
t and f, 


where 
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• for every 0 < k < n — 1, Pk{R, S,T) is homogeneous of degree 2 in the variables 
R, homogeneous of degree n — l — k in the variables S, homogeneous of degree 1 in 
the T variables, 

• in the case k = 0, one has 


Pq {{Rj)j^ {Ta,g)a,g) — ( ^ ^ ^ ^ £{o')RpR(T{p) | 

Vp=l cr€Sn jj^p / 

At first sight, this lemma seems rather technical; yet, writing the integrals this way 
will simplify the presentation afterwards. This lemma gives at least a formal expansion in 
powers of A“^, and we will verify in the following paragraphs that each term in the sum is 
in fact integrable and has uniformly bounded norm in R\. 


Proof. Write hrst that 


1 


lN*i{f=0}) 


u{y,dy,r],dr]) = ^ &\f,f)dt, 


p=i 


where 


dyJ 


&\f, t) = ( Ga(/)*(5”+^), ^df ^|\d[ -dyJ - r/, ) A u{y, r/, dy, dr,) 


j¥=p 


A 


T*U 


Expanding the df term in the wedge product, one can also write 


1 

An 


/ ^{y,dy,ri,dri)= ^ (- 1 )’^+^’ ^ / &''^\f,t)dt, 

W*({/ = 0}) p g^l JM.+ 


where 


t) := / GA(/)*(^r‘). A ( Aat - x) ^ “(S'' 1’ *'> 

Thanks to the dehnition of Ga(/)*(5o^^), this can be rewritten as 

&’‘>\f, t) = ^ (^^^GA{f)*{S^^^),r]pr]g f\d (^-^dyj “ A dyq A u{y, t,, dy, dr,) 


T*U 


Write now 


A%, 

\a\ + \y\=n 


u : = 
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where is compactly supported in T*U. By expanding the wedge product involved in 
the integrals, we can rearrange the sum in the following way: 


1 


A 


' i^{y,dy,r],dr]) 

N*{{f=0}) 
n—1 

k=0 


C/xI 


where 

• for every 0 < fc < n — 1, Pk{R, S, T) is a polynomial which does not depend on A, 
t and /, 

• for every 0 < fc < n — 1, Pk{R, S,T) is homogeneous of degree 2 in the variables R, 
homogeneous of degree n — 1 — A: in the variables S, and homogeneous of degree 1 
in the T variables. 

• in the case A: = 0, one has 

Pq {{Rj)jy (Ta,l3)a,f}) — I ^ ^ ^ ^ ^ip)RpRcr{p) J_ J_ ^j,a{j) 1 • 

\p=l <T&Sn j^P J 


This concludes the proof of the lemma. 


□ 


4.2. Notations and conventions. In the following, we £x the following convention, for 
every 0 < A; < n — 1, 

In particular, one has 


n—1 


(26) J(f,u;,U):= — 


Y‘({/=0}) 


uj{y, dy, r],dr]) = - ^A uj, U). 


k=0 


In the following, the letter J (resp. A', resp. L) will denote functionals of (y, rj, t, f) that 
have been integrated against the variables {y,ri,t) (resp. {y,ri,t,f), resp. /). 


4.3. Integrability. For any element oj := \+\y\=n‘^o^Ay^v)dy'" A dyp in VniT*U), let 

llcull = sup ^ \ua,i3{y,r])\. 

a,y,{y,v)ef*u 

Using the conventions of paragraph 14.21 we will now prove the hrst part of Theorem 14.11 
which follows immediately from the next proposition 

Proposition 4.4. Let u be an element in P„(T*f/). Then, for every 0 < k < n — 1, one 
has 
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belongs to cI^a) for every compact set JC, there exists a constant C/c^k > 0 such 

that for any test form uj in Vn{T*U) supported in 1C, 


'Qa 


\J^’^\f,u,U)\d^IAif)<C^,k\\uJ\ 


Showing this proposition is our more delicate task from the analytical point of view as 
we have to justify carefully several inversions and convergences of integrals. 

Remark 4.5. As was already mentionned, a direct consequence of this lemma is that, for 
every uj in Vn{f*U), 


f eA —— 

An 


/ ^{y,dy,v,dv) 

lN*{{f=0}) 

belongs to L^{'HA,dyiA) which is exactly the content of the hrst part of Theorem 14.11 In 
fact, we have something slightly more precise in the sense that 


JUa \JN*i{f=0}) 

Recall from paragraph 14.21 that 


oj{y,dy,p,dv)]dfrA{f) = -A^ / J^^\f,u,U)dfrA{f) + OiA--^). 


'Qa 


jW(/,a;,[/) := e(a) ( GA{fnSlf+^),r-\y.^p) H 


D2 f 

yi ,y‘^U) •> 


p=la£Sn 


d¥=P 


A2 




T*Ux'R* 


In particular, the second part of Theorem 14.11 is reduced to the computation of the ex¬ 
pectation of J^^l{f,uj, U) with respect to the Gaussian measure djjAif)- This computation 
will be performed in paragraphs 14.41 and 14.51 

We will now prove proposition 14.41 in several steps. At several stages of the proof, we 
will use results on pullback of distributions which are discussed in appendix [Bl 


4.3.1. Proof of proposition \4-4\ ^ First step. As a hrst step, we will give a crude upper 
bound on \J^’^l{f,u),U)\. Let a) be a smooth function compactly supported on /C C T*U 
such that, for every multiindices (a, /3) satisfying |q!| -|- |/9| = n, and for every I < p,q < n, 
one has 

V {y,r]) e f*U, \uJaAy^v)VpVq\ ^ d}{y,r]). 

Moreover, we choose u so that sup^y u{y,7]) < C'||a;|| where G > 0 is a constant 

depending only on the support of u, and where ||.|| is the norm on smooth forms dehned 
at the beginning of paragraph 14.31 As explained in Remark FB.91 Ga(/)*((5o^^) is a positive 
distribution on T*U x R* hence a positive Radon measure GA{f)*{SA^){d^y,d'^r],dt) on 
r*17 X R*. In particular, one has 
(27) 


'T*Ux'R* 
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for some homogeneous polynomial Qk of degree n — 1 — k with nonnegative coefficients 
and which is independent of u, U, {y, rj, t) and A. Then, according to flTTjl and to propo¬ 
sition [SHOl the integrability of U) with respect to dpA will follow from the inte- 

grability of the map 


(y,r,,t,!) e ru X B* X ((|Srt,//A"|) 

with respect to the Radon measure , dl^y, d^rj, dt) where 


, A'(A)ll/ll^ 

hAl/):= e 2 


( N{A) \ 


N(A) 

2 


V 271 y 

Thus, according to (1271) and in order to prove the lemma, it remains to show that 

Jf*UxR*xQA ^ ^ 

is uniformly bounded with respect to A. Let j := (ji,72, ■ ■ ■, j2(n-fc)-3, J2(n-fc-i)) be an 
element in {1,... is an homogeneous polynomial of degree n — 1 — A;, it 

is sufficient to prove that 
(28) 


n—l—k 


I w(9,i,)|ir-’-‘ IT 

r*C/xK*xOA p=i 


K(^Hj,u;,U) := 

is uniformly bounded with respect to A 


aA f 

'-'yJ2p-ly32pJ 


h? 


yiKU)G*M^^){di\d-y,d^r^,dt) 


4.3.2. Proof of proposition If.f] - Step 2: Intertwining integrals. We can now use Proposi¬ 
tion [3110] in order to intertwine the order of integration, and to integrate hrst with respect 
to the / variable. For that purpose, we dehne the following nonnegative Radon measure 
on 

= t^A{f)GA{y,v,ty{6Q+^){df), 

and 


n—l—k 


U 


A,2 


(df)= n 


p=i 


y32p-lyl2p 


fiy) 


h? 


tiA{f)GA{y,p,tn6^+y{df). 


Then, combining (|28|) and the “Fubini identity” fl25|) of Proposition I3T0] to Holder inequal¬ 
ity, one hnds that : 




lT*U xR*xOa 


dj{y,v)\t\ 


n—l—k 

n—S—k 

p=l 


d‘^ ■ ■ f 

y32p-l y32p 3 


L 


d}iy,y)\t\ 


n—3—k 


T*UxR* 


P n—l—k 

4 a 


A2 
dA 

y32p—l y32p 


A2 


yA{f)GUSo^y{df,d^y,d^y,dt) 

^ ■■ ) 


yAif)GA{y,V,t)*iSo^^){df) dl^ydJ^ydt 




< 


/ (h(y,r/)|tr-3-y4"f’A (n^)d-yd^ydt. 

Jf*UxR* 


























CONORMAL CYCLE OF RANDOM NODAL SETS 


23 


Finally, we get that : 


(29) (^a)( fiA)(iW^- 


>T*Ux'R' 


The problem is that we do not know a priori if {y, rj, t) i—)■ (hlA) and {y, t], t) i—)■ 

(^a) are integrable functions of {y,ri,t). In particular, the r.h.s. of (l29|) can a 

priori be inhnite. We will now compute (IIa) for g = 1,2 and show that it is in 

fact uniformly bounded in terms of A by an integrable function. We note that the case 
g = 1 is a particular case of the case g = 2 (when k = n — 1). Thus, we only need to 
compute (I^a)- For that purpose, we introduce the following regularization of the 

Dirac distribution, for every hxed f in R*, 


5lt\Fo, F^,..., F^) := 


kod+f^iT 


-g-2m(ro,r).F . ^ 


|n+l 


Remark 4.6. We denote by T = {0} x (R^^^)* the wave front of Using the conventions 
of paragraph 8.2 in |43], belongs to the space ©[.(R^^^). One can verify that (5”t^)m>i 

’ m 

sequentially converges to for the topology on 'Dp(R”’''^) for every hxed f in R*. Com¬ 
bining the sequential continuity of the pull-back operation (Th. 8.2.4 in [43], see also [T7] ) 
to lemma EH we hnd that (Ga(i/, ?7, t)*(4”V))m>i weakly converges to Ga(|/, t)*(5o’'’^) 

’ m 

for the topology on D'(R^(^1) and for every hxed (g, rj, t) in T*U x R*. 


We can then dehne the following sequences of regularized positive Radon measures, for 
every m > 1: 

2 


n—l—k 


(30) 


p=i 


0^ f(v'] 

^yj2p-iyJ2p J yy ) 


A2 


t‘A{f)GA{v,r,,tr{Slf){f)df, 


We will now verify that 


is uniformly bounded in terms of m as follows: 


(31) 


Vf > 0, u 


A,2,m 


m 



where the constant Cq is uniform for (g, -g) in the support of u, for 0 < m < 1 and A > 0 
large enough. Then, combining this to remark [4.61 one can deduce that //^^^^’^^(Da) is 
bounded by the same quantity and thanks to fl29l) . we can conclude that K^’^\},uj,U) is 
uniformly bounded in terms of A which was exactly what we were aiming for. 


4.3.3. Proof of proposition [THI ^ Final step: Gaussian integral bounds on . It 

remains to prove (13T|) . Observe by the dehning equation 


of 4^2!™ fhat 


^A,2,m 


(Oa) 



n—l—k 


n 


df 

y 


,32p-l y32p 


fiy) 


A2 


2 


I^A(f)GA(y,ihty(SpAKdf) 

' m 
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Jn^+n + l 


^n—l—k 

n 

p=i 


2 

j2p-l,j2p 


{G^{y,ri,t).p„){F)(SlV){F)dF 


where we recognize the pushforward measure v\{F) = (Ga(|/, ?7, t)*/XA) {F) on of 

the Gaussian measure /ia described in Proposition [3]9l This allows us to bound AVm (^a) : 


(32) 
where 

(33) 


z/ 


A,2,m 


(Ga) < 


r)drodr, 


Jn + l 


n—l—k 


<^f(^o,^) : = 


n lb -J.- l\j 

/ 2+ , 1 n ^np-id2p^ 

jRn^+n+l 


^-2iF{To,tT).{Fo,...Fn) ~y 




and f\{F) is the Gaussian measure of Proposition 13.91 This quantity can also be rewritten 
more explicitely using Wick’s lemma as 


(34) 5f(ro,r) = (-l) 


n—l—k 


n—l—k 

n 

p=i 


d 


2FdT, 


J2p-1 J2p 


2-^^(T,A^(y)T) 

Volg(M) 


|T=(To,iT,0) 


where A\{y) was dehned in proposition 13.91 Equivalently, one has 

for some polynom Ry ^ depending only on y, A and j and whose coefficients are uniformly 

bounded in terms of A and y. Gombining this to fl32|) . one finds that 

(35) 


F 


A,2,m 


(Ga) < 


-M 


1 d 


t^^^detBA{y) \2m dy'^ \2m dy[ 


1 d 


KKn, 


^-w{y',BA,miy)v') 


\y'={o,v/t) 


where BA,m{y) = vo 1 ^(m) ^a L) + An' From this expression, we deduce that there exists 

a constant Co > 0, for every {y, rj) in the support of u (which is a compact subset of T*U) 
and for every m > 1, one has 

Vf > 0, 4Vm(^A) < 

which is exactly the upper bound fl3T]) from which we get : 


/r*f/xR* 

^ Ga'IIuZ 


|tp-3-fcle 


<+oo 


for some constant Ck: depending only on the compact support /C C T*U of oj. The 
conclusion of the proposition 14.41 follows. 
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4.3.4. Consequence of the integrahility. We fix /c = 0. We note that we have in fact proven 
something slightly stronger than what we claim in the proposition. Namely, the map 


{y, V, t, f) ^ r \r]„(p) 


cP- f' 

A2 




belongs to L^{T*U x M* x thanks to Proposition 13.101 we can integrate 

this quantity in any order, i.e. either hrst w.r.t. / (and then [y, rj, t)) or hrst w.r.t. {y, rj, t) 
(and then /). 


4.4. Expectation of the leading term. In the previous paragraph, we have verihed that 
all the terms in the asymptotic expansion 




1 

An 


W‘({/=0}) 


u{y,dy,r],dr]) 


n—1 

k=0 


belong to d^pf) and that their norm is uniformly bounded with respect to A. We 

will now compute precisely the expectation of the leading term J^^\f,uj,U) from which 
the second part of Theorem 14.11 will follow. 

Recall from remarks 14.51 and paragraph 14.3.41 that the expectation of the leading term 
decomposes as 

/ n 

EE e{a)Kp^^{uj,U), 

p=l aeSn 


where 


P^p,cr (x^) U'^ 


IT*UxR* 


^hpha(p)^(0,...,0),(i,...,i) {y, v)Lp,a{y, V, t)dJ^ydJ^ 7 ]dt, 


and 

LpAy^v^i) ■= yAf)GK{y,ycy{^A^){df)- 

In particular, it remains to compute the value of Kp^„{uj, U) and before that the value of 

LpAy^vy)- 


Remark 4.7. If one wants to compute a Gaussian integral like 

[ ( [ C?hA(/), 

JHa \JN*{{f=0}) J 

then a natural approach is to use the Kac-Rice formula which follows from the coarea 
formula - see e.g. Th. 4.2 in [TT] or Appendix C in [50]. This kind of formula allows to 
“intertwine” the order of integration in our integral and to integrate first with respect to 
the Gaussian variable. As far as we know, there is not (at least explicitely) a Kac-Rice 
formula for the conormal cycle avalaible in the literature even if this kind of approach 
should a priori allow one (modulo some work) to “intertwine” the order of integration. 
Here, we decided to use an alternative approach to treat this issue which is more based on 












26 


NGUYEN VIET DANG AND GABRIEL RIVIERE 


microlocal techniques such as the pull-back theorem for distributions - see e.g. Appendix [Bl 
This point of view of course also leads to an “inversion” of the order of integration as shown 
by 


4.4.1. Regularization of Kp^u{uj). As in the proof of Proposition 14.41 we will regularize the 
distributions in order to proceed to the computation of U). In the notations 

of paragraph 14.3.31 we have : 


Proposition 4.8. Consider the sequence of approximations of . For every a in 

’ m 

Sn, 1 <p <n and for every fixed {y, p, t) in T*U x R*, set 


n+l 


Then (?/, p, t) equals the integral expression : 

1 d 

.'e ^ I I ( ; 

Jn+l 


n 


2TTdF 


j.o-(i) 


27y^{T,AAiy)T) \ 
g Volg 


voig(M) j dr^dr 

J |T=(ro,tT,0) 


Proof. By the argument of paragraph 14.3.31 we know that 

^n—l—k 

n Llu) I i^l(F)(Sr^')(F)dF. 


AA(/)n%T^G'A(»,',,*)•(«')(#)= 


A2 


Jn^+n+l 


P=1 


Ha j^p 

Since h\{F) is a Gaussian measure with covariance Aa, application of the Wick Lemma 
yields that /ir„ 2 +„+i (UpZl~’' Fla{j)) hl{F){5^^^){F)dF equals 


e 'g m 


in + l 


n 

jAp 


d 


2% dTj^cj(^j) 


27r2(T,AAfa)T) 

Vol£,(M) 


drodr 


\T={TO,tT,0) 


which concludes the proof. 
Note that 


□ 


/-(n 

VjAp 


a2 f' 

yPy'^U) J 


A2 


Rk{f) 


tends to 0 as ||/|| tends to inhnity; thus, it can be approximated (in the topology) by 
a function '0i(/) in C“('Ha)- Moreover, mimicking the final step of the proof of proposi¬ 
tion 03] (namely paragraph 14. 3. 3p . one can verify that the sequence 


'Ha 


yMrGdy,g,ty(pV){df) 


m>l 


is uniformly bounded in terms of m > 1 (not necessarily in terms of A, (?/, p,t)). Combining 
these two observations, we find that, for every 1 < p < n, for every a in Sn, and for a.e. 
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(?/,? 7 , t) in T*U X M*, 

r d‘^ f 

where the last equality follows from fl60|) from the appendix. Mimicking one more time the 
proof of paragraph 14.3.31 one also knows that there exists a constant Co > 0 such that, for 
every m > 1, for every 1 < p < n, for every cr in S'„, one has 

V(?/,r/,f) e supp(a;(o,...,o),(i,...,i)) x M*, r/, f)| < ^e"^. 

In particular, from the dominated convergence theorem, we can deduce that 


(37) 


K, 


P,<7 


OJ 


= lim 

m—>-+00 


j.n —3 


'T*UxR* 


VpVa(p)UJ{0,...,0),{i,...,i){y, v)L^r!S(y^ t)d^yd^vdt. 


4.4.2. Computing the limit asm ^ +oo. Using Gaussian integration, we will now compute 
explicitely the limit appearing in (l37|) as m tends to +oo and as A —)■ +cx). According to 
Proposition 14.81 one has 




g 'g m 

Jn+1 


n 


d 


271 dT, 


j,^U) 


2 n^(T,AA(y)T) \ 


j dr^dr. 

J |T=(T0,iT,0) 


Recall that the expression of Aa(p) was given in Propostion l3.91 If we plug in this expression 
into the previous equality, one can also verify that there exists a constant Co > 0, such 
that for every (p, rj) in the support of a) and for every m > 1, one has 


(38) 

where 


W>0, 


C„A-‘ 


P,(J 




(- 1 ) 


n —1 
2 


(27i)'^ ^ J^n+1 


e e ™ 


n 

j¥=p 


d 


-e ^ 0 * 


{T,AqT) \ 


Vol, 


,(M) 


/ T= 


drodr. 


T=(To,tr,0) 


In particular, provided that the limit limm^+oo {r],t) exists a.e., one can combine (|37|l 

to the dominated convergence theorem in order to get 

(39) 




\Va{p)UJ{ 0 ,..., 0 ),{i,...,i)iy,v) ( lim d^yd""vdt + 0{A ^). 

\ m—^+oo / 


JT*UxR* 

Remark 4.9. We note that, thanks to the particular form of the matrix Aq, one has 

(-1)^ /• ,,.,2. *2,,,2 


where 


iS (>).*) = 


Hp.aiCo) ■ 


^2l'n:T.p- ' ^ Vol„('MUr., 


(21,] 

n 


n—1 


Jn + l 


g Volg(M)(n+2) 


d 


dT, 




-^7^(ToMTi,j)),AoiToMTi,j)) 


r=(To;0;0) 
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After integrating over the r variable and taking the limit m —)■ +oo, one hnds 

n lx \ ^ 

(- 1 )~ 1 /(n + 2 )Volg(M)\ 2 (n+2)Volg(M)|„|2 


lim = 


m —>-+00 


(27r)’^-i r 


27r 


Hp^a{To)dTo. 


4.4.3. Summing over p and a. Recall that everything we want to calcnlate is the leading 
term - J^^J^^\f,u,U)dpA{f) of ^ a;(i/, di/, ? 7 , d?]) which, thanks to ([32]) and 

to Remark 14.91 is given by 

/ n 

ro'>u,u,u)d^Af) = Y.Y. 

p=l (xeSn 

i’""^dpdfT(p)C^(0,...,0),(i,...,i)(d, d) ^)) d^yd^Tgdt 


'T*Ux'R* 


p=l (TGSn 

r) — 1 

(-1)— ('(n + 2)Vol5(M)'^^ 

W. 


n—1 


2tt 


Y. Y 

p=l (tG5„ 


'T*r/xR* 


dpdo-(p)^(0,...,0),(l,...,l)(l/5 d) ^3 


(n+2)Volg(M)|7)|2 

2^2 


Hp^^{TQ)dTQd'^yd'^pdt 
JK. 

In order to alleviate notations, we introdnce: 

n n, 

(40) L(r,) = EE ^(o')dpd<7(p) / Hp^airo)dTo. 

P=i cres™ 

Then, we express everything in terms of L{p): 


’Ha 


(-1)^ / (n + 2)Volg(M )y 

(27r)”“^ \ 27r 


p=l (tGS„ 


'T*U xR* 


IV(0,...,0),(1,...,1) (d) d)dpdo-(p) 


1 _(n+2)Volg(M)|,7|2 


Hp^a{To)dTod^yd^pdt 


2ti 


(-1)"^ /(n + 2)Vol3(M) 
(27r)”-i 
1 

(-1)"^ /(n + 2)Vol3(M) 


/T*t/xR* 


lV( 0 ,..., 0 ),(l,...,l)(l/ 5 d) 


1 (n+2)Vol3(M)|7)|2 

—e 2 p L{p)d'^yd'^pdt by dehnition of L{p) 


(27r)”-i 
X 


27r 


+°° (~+2)Volg(M)|g|2t2 dt 


ee- 


f2 


/ C^(0,...,0),(!,...,!) (d, d)^(d)c^”dc(”d 

Jt*u 

C~°° , (n+2)Volg(M)|g|2t2 

+ / 5 


'0 
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(-1)"^ /(n + 2)Volg(M) 


(27r)^-i 

^+oo 

/ 


27r 


lT*U 




(n+2)Volg(M)\7]\^t^ 


2tdt 


2(-l)^ {{n + 2)\o\g{M))^ 


--1 


{2tiY-^ 


(27r): 


I (^{o,...,o),{i,...,i){y,v)^^d'^yd^r]. 


Therefore, 

(41) 

' 1 


'Wa 


A^' 


'n*({/= 0 }) 


uj{y, dy, y, dy) dy^if) = -4 


[, T^dy^A.. .dy'^Au{y,dy,V,dv)+^i^ 

'f*u Ihl 


where 

2(-l)°^((n + 2)Vol,(M))*-‘ 

"■ (27r)“-i ( 2 jr)i 

In order to conclude, it remains to compute L{y). This combinatorial calculation will be 
the purpose of the next two paragraphs. 


4.5. Computation of L{y). Recall that L{y) is given by the sum fl40|) . We start by 
computing explicitely each term in the sum. We £x a in Sn and 1 < p < n. Moreover, we 
set 

which appears in the dehnition of Hp^„{To) - see Remark 14.91 Our hrst step will be to use 
the Faa di Bruno’s formula for partial derivatives |10] to give an explicit expression of 



2tA {f,Anf) 

Volg{M) 


Recall that this formula states that 



27v'^{f,Ar,f) \ 
Volg(M) j 


27r2<f’,AnT> 

Volg(M) 


E 


/ 27r2 aH(T,yloT) 

Vol,(M)y' 


where 


is the set of partitions of {1 < j < n : j 7 ^ p}. We now £x some notations: 


'Pa,p ■= {21 partition of {1 < j < n : j 7^ p} : Vo G 2t, |o| < 2 and cr(o) = a} . 


Remark 4.10. We note that is empty if cr(p) 7 ^ p. Moreover, for every a such that 
a{p) = p, the above subset is empty if a is not a product of disjoint 2 -cycles. 


For a given 21 in Po-,p, we dehne 

• 2li := (a G 21 : |o| = 1}; 

• 2 I 2 := (a G 2t : |o| = 2 and a\a = id}; 

• 2 I 3 := (a G 21 : |o| = 2 and a\a Y id}- 
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Using these conventions and the specihc form of the matrix Aq - see paragraph 14.3.31 one 
hnds that 


n 


d 


dT, 




27r^{T,AQT) 


g Vo\g{M) 


^2^2 

IQ 


_ g Volg(M) 






Volg(M)(n + 2)/ \n + 4 


|T=(ro;0;0) 

|21| / 1 X |2t2| + |a3| 




Observe now that 

r 27r2-r| 

/ Tq ^ e = 




2ir 


dvr^ 


if |2li| = 0 (mod 2), and that it is eqnal to 0 otherwise. We now set 


n ■ — 




2tG'Pa,p:|ai|=0(mod 2) 


n + 2 


|2l| 


/ 4 2 \ | 2 l|-^ / ^ X 

(^) 


X |2i2|+|a3| 


whenever a{p) = p and C'o-,p := 0 otherwise. 

Remark d.ll. We note that the snm is empty as soon as n is even. In fact, reqniring 
|2li| = 0(mod 2) imposes that n — 1 is even. 


From the above calcnlation, we deduce 

n „ 

(42) L(„) = 5^5; / Hp^aiTo)dTo 

P=i (Te5„, 


Vol,(M) 


27r 


p=l cr^Sn 


In particular, from remark Id. Ill one has L{p) = 0 if n is even. We also note that Xlo-gs ^<^,p 
is independent of the choice of p. Thus, one hnds, whenever n is odd. 


(43) m = 
From (ITTil . one hnally obtains 

(44) f 


27r 


1^1^ Y1 

a&Sn 


'Ha 


1 

A" 


/ (^{y,V,dy,dr])] dpAif) 

/a*({/=o}) / 


Volg(M) 

27r 


^n( / dy^ 

V-.e5„ J df*u 


A.-.dy"^ Auj{y,dy,p,dp) + 0{A ), 


where An was dehned in fldTj) . 
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4.6. Computation of fj44jl . it remains to compute ''^hich 

is independent of the choice of p. For that purpose, we £x 0 < /c < and we suppose 
that a is the product of k disjoint 2-cycles cxi,..., < 7 ^. Then, we compute the value of 
in this case. We note that in this case [Slsl is necessarily equal to fc, and thus 


r — 

^cr,p 


n -I- 4 




u -|- 2 


(|a,|-i)!! 


47r 


2 \ | 2 I|-!¥ 


Vol,(M) 


j 


aeRa,p:|aiNo(mod 2) 

We also note the following useful relations 

fi — 1 = 1 1 T 21 ^2 1 T 21 ^ 3 1 and = \ \ \ ^[2 \ + | *^ 3 1. 


u -|- 4 


I 2 I 2 I 


Recall now that the number of permutations in {1, 2,..., 21} with I cycles of length 2 is 
equal to - see Th. 6.9 in [15]. From this observation, we deduce that 


C = ( 

^ (n + 2)2Volg(M) 


n —1 
2 


n -I- 2 
n -|- 4 


n —1 

k 2 


-k 




1=0 


n + 2 
n -|- 4 


and thus 


27r2 


\2(n + 2)2Volg(M) 


n —1 
2 


2 (u-f 2 ) 


n —1 

k 2 


—k 


E 


(n — 1 — 2k)\ 


^ + 4 J ^ l\ - {k + l))\ \ n + 4j 


n + 2 


Thanks to our formula for Ca,p when a is a product of k disjoint 2-cycles (thus e{a) = 
(— 1 )^), we deduce that 


£{(^)Ca,p 

cr&Sn 


27r2 


n —1 
2 


n—1 n—1 

2 /0 7.M 2 

'i2k 

{n + 2y^Vol,(M) J 


-k 


Ec'f-.w2‘ E (-1) 


I (n — 1 — 2k)\ f n + 2 


l\2^ - {k + l))\ \n + 4: 


l+k 


After simplification, it gives 


Y = in - 1)! 

O-eSn 

and thus 


27r2 


(n + 2)2Volg(M) 


_ , n-l p 

2 / I o \ ^ 2 

-I- 2 ' 


E 

k=0 


E 


(_n±2\ 

V n+i) 


^ + 4 ; ^ kW. - {k + l))V 


Y ^{(^)Ca,p = (n- 1 )! 


C^Sn 


27r^ 


E 


(u-F 2)2Volg(M)/ ^ k\ - k)\ \{n + 4) J \ u-1-4 


n -I- 2 


1 - 


u-F 2 \ 2 




Finally, we obtain 

(45) 5^ E{a)Cy, 

cr£S„ 


[n — 1)! 

TW 


27r^ 


(n + 2)2Vol3(M) 


n—1 

2 
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4.7. The conclusion. Regarding fHB . fl44|) and 

2 (-l)^(n-l)! TT^ 


Bn : = 


, we define 

/ T X n+i 
(- 1 ) 2 


(27r)—1 (n^)! 2(n + 2)t 7rVol(§—i) (n + 2)f 

if n = 1 (mod 2) and = 0 otherwise. Finally, still thanks to fl44|) . we conclnde 
1 f 


I'Hk 




'y*({/=o}) 


uj{y, dy, T], dy) dyA{f) = Bn v)d'^yd'^r] + 0{A ), 


'T-U 


which is exactly the content of the second part of Theorem 14.11 


5. Formal derivation of the main result via Berezin integration 

The derivation of the main theorem (namely Theo. 14.11) is slightly technical dne to the 
fact that we have to justify carefully the smallness of several terms and the inversion of 
several integrals but also due to some complicated combinatorics. The purpose of this 
last section is to show how the combinatorial aspects of the proof can be treated more 
easily at the expense of a little bit of abstraction by making use of the so-called Berezin 
integral |5l EH |25l [63] . 

In this section, we will not pay too much attention to the inversion of integrals and 
to the size of the remainder terms, and we will mostly focus on the computation of the 
leading term of the asymptotics. Similar arguments as the ones presented in the previous 
section would in fact provide a rigorous treatment of the calculation developped in the 
present section. The reason for presenting this alternative approach in a formal (but non 
completely rigorous manner) is that we believe that it can be helpful (at least for the 
reader familiar with Berezin integration) to understand and follow the main lines of our 
proof, especially its combinatorial aspects. It is also plausible that this kind of intuitive 
calculation could be used for other related questions. 

5.1. A brief reminder on Berezin integration. Berezin integration is a convenient 
formalism that allows us to represent the wedge product appearing in our formulas as 
oscillatory integrals. This formalism was introduced by Berezin [S] and it is a more or less 
standard tool in quantum field theory where one aims at computing integrals over both 
bosonic (even) and fermionic (odd) variables. Basic introductions to this formalism are 
described by Disertori [25] or by Tao [63] with a view towards random matrix theory, a 
very good introduction can be found in the fantastic book of Takhtajan [HU Chapter 7] - 
see also the book of Guillemin and Sternberg [37] for considerations on Fourier transforms 
in this context. Following closely the presentation of these references, we give here a brief 
(and somewhat simplified) overview on it that should be sufficient for the purpose of our 
formal calculation. 

5.1.1. Odd and even variables. Let R be a finite dimensional real vector space. One can 

define the so called exterior algebra AV := A^V. For simplicity of notations, we will 

denote by $i<F 2 := <hi A<F 2 the product on this algebra. This algebra can be splitted in two 
parts. One part is made of the so-called even (or bosonic) elements ARven := 
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while the other is made of the so-called odd (or fermionic) elements A 14 dd := 0 ^^ 

Any odd element fl commutes with every even element p, while it anticommutes with odd 
elements. From this observation, one can deduce that any element <F in is 

nilpotent. One can verify that any element <F of NV can be exponentiated as follows: 

<|,fc 

exp(«>) 

/c=0 

Important properties of the exponential are that exp(11) = 1 + 11 for any odd element IT 
and that, for any even element p, 

(46) exp(p + 4)) = exp(p) exp(<I)) = exp(<I)) exp(p). 

5.1.2. Integration of functions of odd and even variables. One of the aim of Berezin in¬ 
tegration is to integrate function whose variables are in AV. In the following, we just 
need a simple version of this formalism. Namely, we fix two nonnegative integers di and 
^ 2 , and we consider the case where 1/ = M. We are only aiming at the integration of 
functions depending on the even variables (pi,... ,Pdi) G (A'^K)'^^ ~ and on the odd 
variables (IIi,..., n^j) G Said differently, we consider the algebra C'°°(M'^^) whose 

elements are smooth functions /(pi,... ,Pdi) of the even variables (pi,... ,Pdi) then we con¬ 
sider the Grassmann algebra G°°(M'^^)[ni,..., If^j] which is a polynomial algebra generated 
by (Hi,..., Ildj) satisfying the relations Iljllj + Iljllj = 0, V(A j) G { 1 ,..., ^ 2 }^ see (HU 
Chapter 7]. More precisely, our functions will be of the form 

f(p, n) := /„(p) + E E A.-3?)^,... n... 

fc>l 

where /* are functions defined on and with values in AV. Given such a function /, we 
will say that it is defined on its Berezin integral is defined by the fundamental 

formula : 

(47) [ /(p, Il)dpdll := [ du^ ■ ■ ■ fdpidp2 ... dpd ^. 

We refer the reader to [25] for a more general presentation of the Berezin integral, for basic 
properties on this integration procedure and for a review of simple examples. 

Remark 5.1. Following [37] - Chapter 7, one can also define a notion of “super Fourier 
transform”, as follows: 

Vg G {A^RY\ vg G (A^RY^ Q) ■= [ exp (- 2 * 7 r(g.p + Q.B)) f{p, n)dpdn. 

jR(‘ill'*2) 

The prefix “super” is here to emphasize that we consider functions depending both on even 
(bosonic) and odd (fermionic) variables. We observe that q.p + g.II is an even element. 

Note that any one form Q can be written as follows: 

(48) g = — 77 ^ f exp(— 2 i 7 rng)(in, 

2*71" Jr(o|i) 
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and more generally any product of d 2 one forms Qi,, Qd^ as 

1 


(49) 


Qi ■ ■ ■ Qd2 




[ exp I — 2 z 7 r IljQj ) dll 

jR{0\d2) V ^ y 


1 ... dU 


d2 ■ 


This formula gives us a representation of products of odd elements as “oscillatory” integrals. 
This observation is the first key point which is at the heart of the calculation below. 


5.1.3. Representation of integration currents as oscillatory integrals over even and odd 
variables. Using the conventions of Appendix |Dl we consider a submanifold S in some 
larger manifold X given by equations G = {gi)i^i^d = 0. Then it is a well known fact that 
the delta function Sq o G = G*(Sq) supported on S' can be represented as an oscillatory 
integral. So a natural question is how to represent the integration current [S'] = [{G = 0}] 
as an oscillatory integral ? The answer involves the introduction of odd variables and the 
use of the Berezin integral. 

Proposition 5.2. Let X be a smooth oriented manifold with orientation [cui], 

G ■. X e X ^ y = {gi{x),..., gd{x)) G 

a smooth function such that the differentials dxgt are linearly independent for every x E X 
satisfying gi{x) = ■ ■ ■ = gd{x) = 0. Set S to be the submanifold in X defined by the regular 
system of eguations {gi = ■ • ■ = fi'd = 0} and oriented by [oj-fi. If dgi A • • ■ A dgd has 
orientation compatible with the pair of orientations [cui], [ufl then the integration current 
[S'] is represented by the oscillatory integral 

(50) [^] = G*( . Pdy^+^ddy^)dp,... dpddU,... dBrf) . 

\{-2z7iY J^(d\d) J 

Proof. The proof follows from a simple calculation, we start from the representation formula 
for the current of integration given in appendix |D] : 

[^] = G* {S^dy^ A ... A dy^) . 

Then, using the classical representation of the delta function by oscillatory integral, one 
has 

G*S^G* {dy^ A ... A dy^) = G* (^j dpi... dp^ G* {dy^ A ... A dy^^) 

and, using the representation flT^ of the exterior product dy^ A • • • A dy^^ in terms of the 
Berezin integral we end up with the expected formula. □ 

In the case where the ambient manifold X is (e.g. in a local chart), we can in fact 
write (at least formally) 

[>^] = r4-Xd f ... dpddU, ...dUd 

[—AtTlj J^(d\d) 

This oscillatory integral formula is already known in the literature since it appears in the 
work of Frenkel, Losev and Nekrasov on instantonic quantum field theory [301 P- 23]. 
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5.1.4. Fourier transform of Gaussian integrals. We observe that the following holds: 

Lemma 5.3. Let N > 1 and A be a symmetric, positive definite N x N matrix. Let 
p = {pi,... ,pn) be an element in {NVeven)^■ Then, one has 


(51) 


{p\A Ip) 

e 2 


det(A); 

(27r)T 


,Ar 




where {c,p) = and {p\A ^\p) = Y.i,jPiPj{A ^)ij. 

This lemma will be the second key point in the compntation of the asymptotics via 
Berezin integration. 


Proof. Let C[[Xi,... ,Xjv]] be the commntative algebra of formal power series in N inde- 
terminates (Xi,..., Xat). First, we note that 


(52) 


and 


det(A) 

(2Tr 


jy. 

2 


ilex') —LLdM ,N 


are well defined as formal power series since the integral on the right hand side is defined 
by expanding the exponential function p-AA powers of X under the integral sign where, 
for each (m, ..., 4 ), the term /]giv(cii... Ci^)e 2 d^c is a convergent integral that can be 
explicitely calculated by the Wick lemma. Then the equality between the two formal series 
in (|52il follows from the classical result on the Fourier transform of the Gaussian measure 
replacing the formal indeterminate X by ?/ G and by uniqueness of the Taylor series 
expansion m. y aX y = t). □ 


Remark 5.4. Conceptually, we note that this proof shows that we can define the Fourier 
transform of a Gaussian measure in a purely algebraic language. 


5.2. Asymptotic expansion in a local chart. We fix a local geodesic normal coordinate 
chart (G, 0) and u){y,dy,r],dri) an element in Vn{T*U). Recall that we want to compute 
the expectation of 

■Jif,aJ,U) := / uj{y,dy,r],dri) 

jN*{{f=0}) 

with respect to the Gaussian measure dfi\. We can write formally 

Recall that one of the reason why the calculation of the expectation against the Gauss¬ 
ian measure is possible is the fact that the Dirac distribution can be represented by an 
oscillatory integral over even variables, i.e. 

Then, some of the (combinatorial) difficulties in the proof come from the fact that the 
wedge product has not at first sight the same simple structure. By Proposition 15.21 we see 


Jif,u,U)= f 

JPUxR* 


^ Au;{y,dy,r],d7]). 
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that the advantage of the Berezin formalism described above is that it allows ns to write 
the wedge prodnct as an oscillatory integral over odd variables in Before writing this 

oscillatory representation - see eqnation (15^ below, we proceed to a few simplihcations that 
will make the calcnlation slightly simpler. More precisely, we expand the wedge prodnct, 
and we obtain that 


r,s=l 


fT*UxR* 


VrVs 


di/* A W^{f, t, y, dy, 17 , dr]) A u{y, dy, r], dr]) A dt, 


where 

t, y, dy, r], dr]) := (f, -r]\ f\d 

For every t in M* and every / in "Ha, this dehnes an element in T>'^^^{T*U). The main 
contribntion in this integral will come from the terms involving derivatives of / - this is 
exactly the content of Prop. 14.41 above. Thns, one can set 

^o{f,'t,y,dy,r]) 




and then one hnds 
(53) 

n 

r,s=l 


r ^r]rr]sdy"AWQ (/, t, y, dy, r])Auj{y, dy, r], dr])Adt+0{K 


n—l\ 


’T*UxR* 


In order to get the conclnsion, everything boils down to compnting the expectation of 
^o{fytyyydy,r]) with respect to d]x\ for a given {y,r],t) in T*U x M*. In order to allevi¬ 
ate notations, we will jnst write W^i^f) and re-establish the dependence in {y,r],t) after. 

Using the representation fj4^ of 1-forms via Berezin integrals, one can write the following 
“oscillatory” integral: 

(54) 

W'l(/) = (-2“ (po/ + P- (^ - p) + E ^ dp,dpdn. 


Remark 5.5. The fact that lTo(/) has a simple representation as an oscillatory integral is 
the central observation that makes onr formal calcnlation mnch simpler than the detailed 
combinatorial proof we gave before. Note that the element Tljdy^ is even since it is the 
prodnct of two odd elements and thns, we are in a sitnation where we can apply lemma [T3l 
As was already mentionned, we do not pay too mnch attention to the inversion of integrals 
and to the size of the remainders, and we mostly focns on the conceptnal aspects that 
make onr proof works. 


Recall that we aim at compnting the expectation of Wq (/) with respect to the Ganssian 
measnre d]XA- For that pnrpose, we £x (po,P, U) in and we nse lemma 15731 to 
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write 


exp -2m pof + p. 


I'Ht, 


tdyf 

A 


E 11 dmf) 


j^r,k 


A2 


(55) = exp 


27r2 


-((po,tp, {^jdy )j^r,k),Ao{po,tp, {Ujdy )j^r,k)) (1 + 0{A )), 


Vol,(M) 

where the covariance matrix Aq was dehned in proposition 13.91 We can now expand 


((Po, tp, (Jljdy’')j^r,k), Ao(po, tp, (Jljdy^)j^r,k)) = 


^ + Po - Y1 


n + 2 


n + 2 


j¥^r 


{n + 2){n + 4) 


{U.dy’^Ukdy^ + U.dymkdy’^). 


j,k^r 


Thanks to the anti-commutation rules for odd variables, one hnds that the last term in 
this sum is in fact equal to 0. Recall that we are in fact interested in computing W^lf) 
which is dehned by (l5T)) . In particular, we have to multiply this Gaussian term by 
and then integrate over the variables (poiP^n) in After integrating over the 

even variables (po,p), one hnds that the following equalities hold: 


e^^P-^exp 


( 66 ) / 

Jr("|o) 

and, by lemma 15.31 


277^^2 IpP 


(n + 2)Vol3(M) 


dp = 


V0\g{M){n + 2)'\^ (n+2)Volg{M)\r,\^ 

^ e 


(57) 


27r2 


/r(i|o) 


exp 


Volg(M) 


Po 


n + 2 


Po 


1 1 ^^0 




VolyiM) 

271 


exp 


27r^ 


Volg(M)(n + 2)2 








It now remains to integrate fl57)) with respect to the odd variables (11^)^^^. Regarding 
the dehnition of the Berezin integral, we have to compute that 


d 


d 


d 


ani ■ ■ ■ dUr 


an. 


exp 


277^ 


Volg(M)(n + 2)2 




j¥=r 


One can verify that if n is even, this quantity vanishes, and then the total expectation 
vanishes. In the case where n is odd, one hnds from the multinomial Newton formula that 


d 


d 


d 


27T^ 


ani ■ ■ ■ an. 


an, 


exp 


Volg(M)(n + 2)2 


E 


j¥=r 


277^ 


(n — 1)! 

(^)! VVol5(M)(n + 2)2 


n — 1 
2 


dy^ A ... dy^ A ... dy"^' 
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Combining this to fl57|) . we obtain 

27r2 


(58) 


(l|n-l) 


exp 


VolgiM) 


pI 


n + 2 


Po ^ Ujdy^ dpodU 


Vo\g{M)Y (n- 1 )! 


27r 


n—1 

2 


j¥=r 

2'k‘^ 


i)! \Volg{M){n + 2f 


dy^ A ... dy^ A ... dy"^. 


From (l56i) and (l58|) . we deduce that 

[ WQ{f,t,y,dy,ri)dpAif) 


'Ha 


(— 1)"2 (n — 1 )! TT 


71 — 3 
2 


(n + 2)Vol,(M)|,7| 


^(n+2)Volg(M)e ^ dy^A... dy^A... dy"'{l+0{A~ )). 


(27r )—1 (i^)! 2 (n + 2 )t 

We now use equality (|53ll . Summing over r and s and integrating over t, one hnally gets 
that, if n is odd, 

2 (-l)^ 


I J (/, (U, U)dpA{f) = u4n-ll ( 

Iha 7 rVol(S" 1) \A/n + 2 


dy^A... dy'^Au{y, dy, p, dp)+0{A 


n—l\ 


fT*U 


Appendix A. Derivatives of the Spegtral projector. 

Recall that (ey)i<y<Ar(A) is an orthonormal basis of "Ha made of eigenfunctions of A. We 
dehne the normalized projection kernel: 


N{A) 


CA{y,z) := 






where N{A) = dim T-La- The reason why we can perform some computations with respect 
to the Gaussian measure pA is that we have very precise asymptotic informations on the 
kernel CA{y,z), at least on the diagonal. Precisely, Bin proved the following result [9] 
building on earlier arguments of Hormander 


Theorem A.l. Let {U,(j)) be a sufficiently small geodesic normal coordinate chart. For 
every multiindices a and /3 m Z", the following estimates hold uniformly for y in (f>{U), 
as A ^ +CX 0 , 

where, for a = (3 = t), 

C - 1 


for a = (3 (mod 2 ) and (a, (3) 7 ^ (0, 0), 
Cn,aP := (- 1 )“ 


Vol,(M)’ 

nLi(«i + - 1 )- 


Volg(M) {\a + (3\ + n)... {n + 2)' 
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and 

Cn,a,i 3 ■= 0 otherwise. 

Remark A.2. We used the following conventions: 

(-1)!! := 1 and (2m - 1)!! := (2m - l)(2m - 3)... 3 x 1. 

We gather the values of the derivatives of the spectral projector which are used to prove 
Proposition 13.91 

Corollary A.3. One has 

• for derivatives of order 0, one has 


\ y=z 


Vol,(M) 


0(A-'), 


for every I < j <n, one has 
1 


^2^yiCAiy,z)\y=, = 

for every 1 < j < n, one has 


(n + 2) Volg(M) 

i < n, one has 

^dy,d,,CA{y,z)\y=, = (^^2)Volg(M) 
for every 1 < j ^ k < n, one has 


C(A-^), 


+ C(A-i), 


^A^9l,CA{y, z)\y=, z)\y=, + 2){n + 4) Vol 3 (M) 

• for every 1 < j < n, one has 


+ C(A-^), 


dyjd^jCA{y,z)\y=, = 


A4 (n + 2)(n + 4)Vol3(M) 

for every a ^ (3 {mod 2), one has 

8 ?afCA(!/, 2 )U = o(A-'). 


0(A-‘), 


AI“+/3| 


Appendix B. Proof of Proposition 13.101 

We £x / in Ha = 'Ha\-Da- Recall that Da is a zero measure subset of Ha and that its 
complement is an open. Given any test function '^{f,y,r],t) in D{Q\ x T*U x M*), our 
goal in this paragraph is to make sense of the following Fubini equality : 

= {l,{GAv.V.tnsri'i)nX.u...- 

The fact that the distribution Ga(/)*^o^^ is well defined for / in Ga follows from the 
following theorem of Hormander - Th. 8.2.4 in [43] : 
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Theorem B.l. Let C and ^2 C be two open subsets. For any smooth map 
G : f2i —)■ ^ 2 ; the normal of G is defined as 

Ng := {(G(x),r) e (^2 X : ro^G = O} , 

the pull-back operation extends uniquely to the distributions u G V'{fil 2 ) whose wavefront 
WF{u) does not intersect the normal Ng of G. Moreover, the wavefront set of G*u is 
contained in the set 

G*WF{u) := {{x,To d,,G) : (G(a:), r) G WF{u)} . 

Remark B.2. For the definition of the wavefront set, we refer to |l3l [16]. Recall that, in 
the case of one has 

WF{5^-^^) := {(0,r) : r ^ 0} C 

In the following, we will apply this Theorem in three distinct situations. Recall that, for 
every / in we defined the following map on "Ha x fifiU) x M"" x M*: 

Ga : (/, y , V , t)enAX fi{U) X M" X M* ^ (^f{y); LdyJ - p ,,..., jdyj - 77,^ G M x 

and our idea is to think of "Ha x 4>{U) x M” x M* as the cartesian product of "Ha with 
(j){U) X M” X M* and to think of the global distribution Ga^q"*"^ ^^so as a distribution on "Ha 
depending smoothly on the parameters in 4>{U) x R"' x R* and conversely as a distribution 
on 4>{U) X R” X R* depending smoothly on the parameters in RLa. To be more precise, we 
will use Hormander’s Theorem in the following three cases : 

(1) We fix / in and we consider the partial map Ga(/) : T*U x R* —)■ R”+^. As 
dyf 7 ^ 0 on /“^(O), one can verify that 

Nom, n WF(sr') = 0. 

It means that if we freeze /, the pull-back Ga(/)*^o^^ is a well-defined distribution 
on f*U X R*. 

(2) We consider the map Ga : LLa x T*U x R;^ —)■ R"+^. In order to apply Hormander’s 
result, we would need to verify 

Ng^ n = 0, 

which is slightly less obvious. 

(3) We fix {y,ri,t) in T*U x R*, and we consider the partial map GA{y,V,t) ■ FLa —t 

situation in this case is also slightly different from the first case. 

The second and the third cases are in some sense related to the notion of ampleness 
appearing for instance in [531150] . and they are contained in the following lemma: 

Lemma B.3. There exists Aq > 0 such that, for every A > Aq, 

Ng^ n WF{S^+^) = 0, 

and, for every {y,r},t) in T*U x R*, one has 
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In particular, defined as a distribution on'HAxT*UxM.* and GA{y,ri,t)*{6Q~^^) 

as a distribution on Ha. 

Remark B.4. We note that in the second and in the third cases, the distributions are 
dehned on the whole space Ha, and not only on Qa- In particular, we can consider their 
restriction to the open subset IIa- 

Before continuing our discussion on the properties of these pulled-back distributions, we 
start by giving the proof of the lemma. 


Proof. Suppose that (0, r) belongs to the normal set NGA{y,v,-t) other case works simi¬ 
larly). Then, for every / in Ha, one has 


Y(A) 

rofiy) + ^ = 0 - 


i=i 

pA fA fA 


Consider now a family of functions /q , ff, ..., /,y in Ha, we write 
det I ( f^{y), ..., ^ = ^det {{f^{y), dyifl ^,..., dy^^ff^) 


A 


0<l<n, 


0<l<n 


If we are able to hnd a family (/;^)o</<n in Ha such that the right-hand side of this 
equality does not vanish, we will deduce that r = 0. For that purpose, we note that we 
can £nc|^ a family of smooth functions /o, /i,..., /„ on M such that the right hand side 
does not vanish for every y in U. Unfortunately, these functions do not a priori belong 
to Ha but we can solve this issue since by Lemma 1C.II proved in appendix, hnite linear 
combinations of eigenfunctions are everywhere dense in C°°{M). Said concretely, to every 
0 < / < n corresponds a sequence (/;^)a in Ha which converges to fi in every norm (this 
is just an approximation property like the Stone Weierstrass Theorem but on Riemannian 
manifolds). From this, we deduce the statement of the lemma for A large enough. □ 


We refer the reader to the section o of the appendix where a similar method yields an 
unusual proof of the Whitney embedding Theorem. 

As a direct consequence of the pull-back theorem of Hormander, we can also describe 
the wavefront set of G'a((^o'''^) as follows: 

Corollary B.5. For A > 0 large enough, the wavefront set of is included in the 

following subset ofT*{HA x T*U x M.*): 

WF{GI{6q^^)) C {{yfr]j,t,f;yj,rift,f) such that GA{f,y,V,t) = {0,... ,0) 

Vj = rodyif + Ti-ajyi/ + • • • + rfj = tj, 

1 n ^ t 

f = - ^ Tjdyif, VI < i < A^(A), fi = Tod A ^ ■^d^.id 

i=i j=i 

for (ro,ri,...,r„) ^ (0, ...,0)}. 


^We may have to pick a slightly smaller neighborhood U. 








42 


NGUYEN VIET DANG AND GABRIEL RIVIERE 


Remark B. 6 . This corollary has the following consequence. The set hhF(G^(^g’'"^)) does 
not contain any element of the form {y\ rjj, t, /; 0 , 0 , 0 , /), with / 7 ^ 0 or rjj, t, /; yj, rj^, t, 0 ) 
with [yj, rp R) 7 ^ 0. The first observation is more or less immediate if one uses Theorem lA.ll 
in the case a = /5 = 0. The second one follows from the proof of lemma lESl 

The relevance of such a remark is due to the following general result which can be found 
in [2ni Proposition 1.3 p. 502] : 

Proposition B.7. Let u be a distribution in D'(M"'i x denote by {x,y,p,T]) G 

T* (M”! X the coordinates in cotangent space. IfWF{u) contains no elements of the 
form such that p 7^ 0}, then u{x,y) is smooth in x with value distribution in y 

in the sense that for all test function (p 2 {y), x 1 —> {u{x ,.), (P 2 (-))k '*2 is smooth in x. 

Thus, one finds that there exists / G Ha Ta(/) G V'{f*U x W) such that, for every 
V’l G T>{Ha) and for every ip 2 ^ T>{f*U x R*), 

and 

{y,yt) ^ {GA{y,V,trS^^\f^ih, e C^{f*U x R*). 

Moreover, one has the “following distributional version of the Fubini Theorem” : 

V(V'i,i/>2) e V{Ha) X V{f*U X R*) 

= R2,{GA(y,ri,t)‘{sy),y ), 

Remark B. 8 . Note that Ta(/) is well defined on the whole space Ha while Ga(/)*((5o’'"^) 
was only defined on IIa by Theorem IB. II For a fixed / in IIa, one has in fact Ta(/) = 

Ga(/)*(C‘). 

By density of the functions of the form V’i(/) M'ip 2 {y, rj, t) (where S is the exterior tensor 
product) in V{VtA x T*U x R*) ~ see Ch. IV in [58], given any test function T(/, |/,r 7 ,t) in 
V{VLa X T*U X R*), one has the following Fubini equality : 

(59) (IXTa, 

= (l,(GA(»,',,i)*('5o”+‘).'I'>«b. „ . . 

\ ' n-A/ T*UxW‘ 

Remark B.9. We emphasize that all the distributions considered so far are constructed from 
the Dirac distribution ^g and that they are positive distributions. In particular, according 
to [58] (Ch. I), they can all be identified with positive Radon measures. Therefore, it makes 
sense to test them against nonnegative measurable functions which are not necessarily 
integrable (or integrable function which are not necessarily smooth). This identification is 
used several times in our proof of proposition I4.4I 
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Remark B.IO. We now make a final useful comment. We fix x 'D(T*U x 

M*). According to fl25l) . one has 


/r*;7xi 


' 4 ’ 2 {y,V,t) ( / ) (Tyd^ydt 


’Ha 


Mf) 


'Ha 


y2{y,V,t)TA{d'^y,d'^V,dt) df. 


’T*UxR* 


From the above discussion, one knows that 


/'-t'0i(/) f [ '4!2{y,V,t)TA{d^y,d'^r],dt)) 

\Jf*UxS.* 

is a smooth compactly supported function. As = 'Ha\Da and Da has zero Lebesgue 
measure, one has then 

y2{y,vR) ( [ i>i{f)GA{y,vRy{SQ^^)idf) \ d^yd^ydt 

CxR* \dHA / 



Mf) 


'Oa 


'T*UxR* 


My,V,t)GA{fy{6^+^){d-y,d-y,dt) df. 


= [ ij2{y,vy)(f yi{f)GA{y,vyy{do^^){df)) d'^yd^T^dt, 

Jf*UxR* XJUa / 

where the second equality follow from the fact that 

G'A(/)*(5o”^')(d>,d>,dt)d/ and Ga(i/, r], 

dehne the same Radon measure on the space Da x T*U x R*. Note also that this equality 
holds for any smooth compactly supported function '^ 2 - In particular, one hnds that, for 
a.e. (y, ? 7 , t) in f*U x R*, 

( 60 ) [ Mf)GA{y,r,,tr{S^+^){df)= f i/>i(/)Ga(»,>?,<)’(C')W). 

J I-La ^ 


Appendix C. A functional analytic proof of Whitney embedding Theorem 

In this short section, we would like to prove briefly the approximation property used 
in the proof of Lemma IB.31 and to illustrate it by giving a simple proof of the Whitney 
embedding Theorem. In fact, on a smooth compact manifold M, if we have enough smooth 
functions ei,... ,eAr so that their linear combinations approximates any smooth function 
with enough accuracy, then we will prove that we can embed the manifold M in R^ for N 
sufficiently large by the simple map : 

a; G M i-A (ei(a;),..., ej ^{ x )) G R'^. 
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C.l. Approximation of smooth functions. We denote by H^{M) the Sobolev space of 

s 

functions or distributions of order s i.e. the set of all distributions t such that Agt belongs 
to with norm ||.||h'*(m) = ||(1 ~ -WL^iM)- Note by the spectral Theorem that 

1. 

( +00 \ 2 

where we use the conventions of the introduction. First, we show that eigenfunctions of the 
Laplace operator span a vector space which is everywhere dense in C°°{M). This property 
was used to prove lemma IB.31 


Lemma C.l. Let {M,g) be a smooth compact boundaryless Riemannian manifold. Then, 
finite combination of eigenfunctions (e,), of the Laplace-Beltrami operator Ag spans a dense 
subspace of C°°{M). More precisely, for every norm ||.||c''' ■' 

N 

ll/“^(/)ei)L2ej||cfe —)■ 0. 


Proof. Eigenfunctions (e*)* of the Laplace operator are smooth by elliptic regularity. Thus, 
any hnite combination of them belongs to C°°{M). By the Sobolev embedding Theorem, 
it suffices to show that for any s > 0, 

N 


11 / - X](/, ei)L^ei\\H- 

i=l 


0 . 

N^oo 


For that purpose, we write 

N 

11 / - ei)L2ei\\]j. 

i=l 


(1+-^?)H(/)ei)L2 

i=N+\ 


which tends to 0 as / is smooth. 


□ 


C.2. Whitney embedding Theorem. We will now show how this approximation prop¬ 
erty can be used to prove the following version of Whitney embedding Theorem: 

Theorem C.2. Let {M,g) be a smooth compact Riemannian manifold without boundary. 
Let {ei)i be the sequence of eigenfunctions of the Laplace-Beltrami operator as defined in 
the introduction. Then, there exists Nq such that for any N ^ Nq, the map 

(61) i]g : X e M (ei(a;),..., cn^x)) G 

is an embedding. 

Proof. We first prove that z at is an immersion then we show it is injective. 

Fix some N. To show that zat is an immersion, we proceed as in the proof of lemma [Bl3l 
It suffices to show that for all x G M, di^ix) has maximal rank in other words 

X HA (dei ,... ,deN) 
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has rank dim M = n. Equivalently, it means that near every xq G M, there is a neighbor¬ 
hood U, n sequences )is:is:Ar, • • •, such that for all x E U, the covectors 

N N 

a\dei{x ),..., a^dei{x) 

i=\ i=l 

are linearly independent. In other words, the n-form 

N N 

a\dei) A ■ • • A a^dci) 

i=\ i=l 

never vanishes on U. But this is easy to achieve since for every xq, we can always hnd some 
sufficiently small neighborhood U of xq and functions /i, ...,/„ so that d/i A ■ ■ ■ A dfn{x) ^ 
0,Vx G U. And since the family (cj)* spans a dense subspace of we can always 

hnd Nq large enough so that the vector space spanned by is very close to /i,..., /„ 

in the topology. 

The last part consists in proving that the above map is injective for N large enough. 
By the previous part, we can £x Nq large enough so that is an immersion. Since an 
immersion is locally injective and M is compact, there is some e > 0, so that for all x G M, 
consider the ball B{x,e) centered around x of radius e, the map iNo ■ B{x,e) i—)■ 1^-^° 
is injective. In other words if iNo{x) = iNoiu) for some {x,y) G then x = y or 
d{x,y) > e. We denote by Am the diagonal in M x M. Consider the map ivTo x ijvo ^ 
(M X M) \ Am I—t x which reaches the diagonal A C x exactly at the 

pair of points {x,y) where fails to be injective. First note that x ipfQ)~^{A) = 
{{x,y) such that = iNo{y)id{x,y) ^ e} is a closed subset of M x M \ Am- Since 

M X M \ Am is precompact, the closed subset {ip^^ x iisfQ)~^{A) is compact. Note by 
definition of iN that for all N, {iN+i x iAr+i)“^(A) C {iN x ii\f)~^{A). Hence consider the 
sequence {{iN x iAr)“^(A))^ for N ^ Nq, it is a decreasing sequence of compact subsets in 
M X M \ {{x,y) G s.t. d{x,y) < e}. By contradiction, assume that for all N ^ Nq, 
the map In fails to be injective, then it means that {In ^ *a)~^(A) is non empty for all N, 
then for all N there is a pair {xn, yn) E M x M \ {{x,y) E s.t. d{x, y) < e} such that 
iN{xN) = iN{yN) for all N. By compactness oi M x M \ {{x,y) E s.t. d{x,y) < e}, 
we may assume that the sequence (xat, yN) n converges to some element {x, y) m. M x M\ 
{{x,y) G s.t. d{x,y) < e}. It follows by definition of iN that ei{xk) = ei{yk) for all 
k ^ i and by continuity of the the functions (cj)* and the fact that Xk ^ x,yk ^ y, we 
also have that ei{x) = ei{y) for all i. 

On the other hand, we know that there is some smooth function / such that f{x) ^ f{y)- 
And since / = ei)ei, this implies that/(x) = 'Et^{f,ei)ei{x) = 'Et^{f,ei)ei{y) = 

f{y) which contradicts the fact that /(x) ^ f{y). □ 

Appendix D. Background material on the theory of currents 

In this appendix, we give formulas for the current of integration over submanifolds 
defined by systems of equations of the form gi = g 2 = ... = Qd = ^ ioi some sufficiently 
nice functions gi, ■ ■ ■ ,gd- This kind of formulas can be found for instance in Schwartz’s 
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book [58] - section V.5 (example 3). The advantage of representing currents of integration 
with Dirac distributions is that they will then be easily integrated against the Gaussian 
measure on Ha - see section 01 For an introduction to the theory of currents, we also 
refer to the books of de Rham [23| and Giaquinta-Modica-Soucek [SB] . 

D.l. Currents on zero sets. In mathematical physics, one is often interested in studying 
subsets in X given by equations {G = 0} for some G = In our article, we will 

encounter two classes of currents associated to this type of submanifolds. The hrst class are 
currents representable by integration on some smooth submanifolds which form a subclass 
of the De Rham chains - see chapter 2 in [21]. Let X be a smooth oriented manifold of 
dimension m, let S' be a smooth, oriented, compact submanifold of X of dimension d < m. 
This defines an integration current [S'] in 'D'^{X) which acts on test forms ca G Vd{X) by 

(62) [S]{u) = [ u. 

Js 

The second class consists in Dirac distributions which are defined as follows. We denote 
by (5g the usual Dirac distribution at {0} on M'^. This is dehned by 

(63) \/ip{x)dx^ A • • • A dx'^ G 'Dd(IR'^), h^^x^ A • ■ • A dx‘^)^d = (p(0). 

Now let us dehne the Dirac distribution G'*((5o) = 6q o G supported by the zeros of a 
function G : X —)■ 

Proposition D.l. Let G = £ G°°{X) be such that, for all x ^ X satisfying 

gi{x) = ■■■ = gd{x) = 0, the linear forms d^gi^x),... ,dxgd{x) are linearly independent. 
Then 

(1) the pull-back of Sq G P(;(R'^) by G which is denoted by G*6 q (or also 6 qO G) is 
well-defined, 

(2) the wave front set of G*Sq is contained in 

N*{{G = 0}) = {(x,0 e T*X : G{x) = 0,e = r o4G for some r ^ 0 G R"^} 

(3) for every seguence (Tp)p which converges to Sq in (R'^), G*Tp G*Sq in 

the normal topology o/P)y,|g,^g|(X) which implies convergence in the weak topology 
ofV\X). 

Remark D.2. For the notion of wave front set and for the precise definition of the set 'Dp(X) 
of distributions whose wave front set is contained in F, we refer the reader to section 8.2 
of [l3] (see also [16]) and for the definition of the normal topology on Pp(X), we refer to 

na. 

Proof. Let Nq = {(?/;h) ^ T*^. s.t. G{x) = y and rj o d^G = 0 for some x G X} be the 
normal set of the map G. By assumption, Nq does not meet the set {0} x (R'^)* C T*R‘^ 
which is the wave front set of 6q. Therefore by the pull-back Theorem of Hormander - 
Th. IB. II below. G*dQ = S'^ o G is well-defined and its wave front set is included G*({0} x 
(R^)*) = Xq({G = 0}). Finally, the last claim follows from the sequential continuity of 
the pull-back operation - Th. 8.2.4 in [43] or [T7]. □ 
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Now the natural question which arises is how the integration current on S' = {G = 0} 
and the delta form G*6q = 6qoG are related to each other. Note that the Dirac distribution 
of the submanifold S = {G = 0} is obviously not the same as the integration current [S'] 
because we integrate top forms against Sq o G but only (m — d) forms on [S']. 

D.2. Orientations on submanifolds. Before comparing these currents, we need to dis¬ 
cuss briefly orientability questions. For any smooth manifold X (orientable or not), there 
exists a smooth principal Z 2 bundle Ox ^ X called the orientation bundle. Recall that 
if X is an orientable manifold, then there exists a volume form u G r(X, A’^T*X) where 
n = dimX or equivalently there exists a global section of Ox- Indeed, every choice of 
volume form u on X defines canonically a global section [ca] of the orientation bundle Ox, 
also [(puj] = [ca] for nW ip > 0. Now, if we are given a submanifold S <Z X dehned by 
equations Qi = ■ ■ ■ = Qd = ^ and if X is oriented by [wi] G Ox and S by [UJ 2 ] G O5, we 
choose representatives oji G r(X, A"T*X) and 0 J 2 £ r(A, A"'“'^T*S'). Then we say that 
the orientation of dgi A • • • A dgd is compatible with the pair of orientations [wi], [uj 2 \ if 
dgi A • • • A dgd A a; 2|5 = puj\\s for some smooth function p s.t. (^[5 > 0. 


D.3. Relation between the different types of currents. In the next proposition, we 
give a concrete description of the action of Dirac distributions G*5q on differential forms. 

Proposition D.3. Let G = {gi)i<^i<^d be a smooth map from X to such that the differ¬ 
entials dxgi are linearly independent for every x satisfying gi{x) = ■ ■ ■ = gd{x) = 0 . 

Then, for any test form u G Vm{X) and for any form a in Vm-d{,X) such that dgi A 
■ ■ ■ A dgd Aa = u, one has 

(64) {G*6iuj)= [ a. 

J { G = 0 } 


Proof. We first prove the claim in the case where g = {gi)i^i^d is a linear map from M” 
to of maximal rank. Then one can choose [gi,... ,gd) to be the hrst d components of 
some system of linear coordinate functions {gi,... ,gd,..., gn) on R”, in this case G is just 
the linear projection R*^ 1 —)■ R"^ on the hrst d components and the zero set {G = 0} is just 
{0} X R”“'^ c R'^. For any function if on R'^, the pull-back of '0 by G simply reads 

G*f’{gu ■■■,9n) = ' 0 (^ 1 , ...,gd)^ Hdd+U ■■■,9n) 

and therefore, by continuity of the pull-back operation [16], one has (5 q o G = (5 q ® 1. Any 
form oj has a unique representation in the basis dgi A • • • A dgn as 

= (^i...n{9i, • • •, 9n)dgi A--- Adgn 


therefore : 


[^0 ° G]{uj) = [S^ ® l(^d+l, . . . ,^n)](^) 

/ ^l...n(07 * * * ; O 5 • • • ; A * * * A 

Note that any form a such that dgi A ■ ■ ■ A dgd Aa = u will coincide with 


a^i...n(0,..., 0, gd-\-iy • • •, gn)dgd-\-i A • • • A dg^i 
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when restricted to {G = 0} = {0} x It follows that the pairing does not depend on 

the choice of a. 

Let us go back to the manifold case. For any m G {G = 0}, let U be some neighbor¬ 
hood of m on which the differentials dgi are linearly independent. Then one can choose 
{gi,...,gd) to be the hrst d components of some system of local coordinate functions 
{gi, ■ ■ ■, gd, ■ ■ ■, gn) on U. What we just did was to dehne some local chart $ : m G 17 C 
X I—)■ {gi{m), ..., gn{m)) G M" on 17 C X such that G o : M" ha is just the linear 
projection on the hrst d components. Then the distribution 6q o (G o <I)“^) is well dehned 
by our previous discussion since G o is linear. Then we have the identity : 

((Got-YiS,u,) = 

Va s.t. dgi o A • • • A dg^ o A a = cj. 

We just constructed a delta distribution (G o <|>“^)* (5q supported by some linear sub¬ 
space and hnally, to go back to the open subset U in the manifold X, we must pull¬ 
back the distribution (G o Sq by $. By the pull-back Theorem IB.11 this yields 
(G o (5q = G*(5q and, for every a s.t. dgi A ■ ■ ■ A dga Aa = u, 

{G*6iu) = 

□ 

In fact, using the notations of paragraph ID.21 concerning the orientation of submanifolds 
dehned by systems of equations, the following corollary holds : 

Corollary D.4. Let X be a smooth oriented manifold with orientation [cji], let 

^ •= : X ^ 

he a smooth function such that the differentials d^gi are linearly independent for every 
X E X satisfying gi{x) = ■ ■ ■ = gd{x) = 0. Set S to be the submanifold in X defined by the 
regular system of eguations {gi = ■ • • = fi'd = 0} and oriented by [ufl- If dgi A ■ ■ ■ A dgd has 
orientation compatible with the pair of orientations [wi], [ufl then one has 

(65) [5] =G*(5o") dgiA---Adgd. 

Proof. We note that this corollary is proved in (Section V.5) when 7=1. It is in 
fact a direct consequence of the property of the Dirac distribution 6q o G established in 
proposition ID.31 i.e. 

A ■ • • A dgd, a) = {G*6q, dgi A ■ ■ ■ A dgd A a) = / a. 

Js 

□ 


/ 

J { G = 0 } 


j 

7{Go<i>-i=o} 
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